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ABSTRACT 

Topological defects are produced during phase transitions in the very early 
Universe. They arise in most unified theories of strong, weak and electromagnetic 
interactions. These lectures focus on the role of topological defects in cosmology, 
with particular emphasis on the models of structure formation based on defects. 
The role of topological defects in baryogenesis is also reviewed. 



★ Based on a Troisieme Cycle lecture series given at the EPFL, Lausanne, Switzerland, 
4/22-5/13 1993. 



1. Introduction 



Topological defects 1 ^ are inevitably produced during phase transitions in the 
very early Universe. One or several types of defects arise in most unified particle 
physics models of strong, weak and electromagnetic interactions. Unless they are 
inflated away 2 -* (i.e. unless there is a period of inflation after the phase transition 
during which the defects are produced), topological defects will play an important 
role in cosmology (for recent reviews see e.g. Refs. 3-6). Main aim of these lectures 
is to give a comprehensive pedagogical review of the role of topological defects in 
cosmology. 

The reason why topological defects can play a role in structure formation in 
the early Universe is very simple. Topological defects carry energy. This energy 
leads to an extra attractive gravitational force, and hence the defects can act as 
seeds for cosmic structures. In particular, cosmic strings (for a recent review see 
e.g. Ref.7) and global textures can lead to attractive scenarios for the formation 
of galaxies and large-scale structure. 

In cosmology there are many open questions. What is the origin of the ob- 
served large-scale structure in the Universe ? 9 ) Why are galaxies clustered? What 
explains the specific amplitudes and slopes of the galaxy and galaxy cluster correla- 
tion functions? 10 -* Do the models of structure formation predict cosmic microwave 
anisotropies in agreement with the recent COBE results 11 )? There is a wealth of 
quantitative data in need of a consistent theoretical explanation. 

Before 1980 there was no causal explanation starting from first principles for 
any of the above questions. Around 1980 there was a crucial realization which has 
now led to the development of several self-consistent theories. The key point was 
to realize that in the very early Universe, matter cannot be described in terms of 
an ideal gas. At very high temperatures and energies, the matter content of the 
Universe must be described in terms of quantum fields. In 'old cosmology' (prior 
to 1980), space-time was described by general relativity and matter as an ideal gas, 
in 'new cosmology' space-time is still described by general relativity, but matter 
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in terms of fields. Evidently, from a conceptual point of view the present theories 
are not satisfactory since they describe one sector in terms of classical physics, the 
other quantum mechanically. 

Since in 'new cosmology' matter is described in terms of fields, the possibility 
of having phase transition in the early Universe arises (for current particle physics 
models of matter this is a rather generic prediction). Hence, topological defects 
are formed, and they may play a role in structure formation. 

Concrete issues which must be addressed and which are discussed in these 
lectures are: 

- Why do defects form? 

- What kind of defects arise? 

- How to they evolve after formation? 

- What are the mechanisms by which defects give rise to cosmic structures? 

- Are there direct observational signatures for topological defect models? 

Phase transition and therefore topological defects arise in many areas of physics. 
In metals and other solids, defects arise during the fluid to solid phase transition. 12 ) 
The resulting crystal defects can be point, line or planar defects. In liquid crys- 
tals, defects arise during the transition from the disordered phase yo the ordered 
phase. 13 ) There have recently been some very ingenious studies of the dynamics 
of defect formation in liquid crystals. 14 ) Vortex lines also appear in the transi- 
tion from the normal to the superfluid phase of the 3 He 15 \ and in the normal to 
superconducting transition in superconducting materials 16 ). 

The above are all examples of topological defects in condensed matter sys- 
tems as opposed to relativistic field theory. There are important similarities and 
differences between the two classes of systems. In condensed matter systems the 
dynamics is generally friction dominated and nonrelativistic, whereas in flat space- 
time the dynamics of topological defects in field theory and cosmology is relativis- 
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tic. Hence, condensed matter systems give a good analogy for the topology and 
formation of defects, but not for the dynamics. 

These lecture notes begin with a survey of the successes and problems of stan- 
dard 'big bang' cosmology. In Chapter 3, phase transitions and the origin and 
nature of topological defects are reviewed. In Chapter 4, the formation and evo- 
lution of topological defects (in particular cosmic strings) in the early Universe 
are treated. Chapters 5-8 deal with topological defect models for the formation of 
structure. First, an introduction to the basic mechanisms and a survey of the most 
relevant data is given. Then, the cosmic string and global texture models of struc- 
ture formation are developed, and in Chapter 8 the predictions of the two models 
are compared with those of inflationary Universe scenarios. Chapters 9 and 10 
deal with the microphysics of cosmic strings. Baryon number violating scattering 
processes based on strings are discussed in Chapter 9, and the role of topological 
defects in baryogenesis is analyzed in the final chapter. These lectures are based in 
part on lectures 1T ) given at the 7 th Swieca Summer School in Particles and Fields, 
(see also Ref. 18). 

A word concerning notation. Unless mentioned otherwise, units in which c = 
h — k,B — 1 are employed. The metric is taken to have signature (+, — , — , — ). 
Greek indices ran over space and time, Latin ones over spatial indices only. The 
Hubble expansion rate is H(t) = a{t)/a{t), where a(t) is the scale factor of a 
Friedmann-Robertson- Walker (FRW) Universe. The present value of H is h x 
lOOkms" 1 Mpc^ 1 , where 0.4 < h < 1. Unless stated otherwise, the value of h is 
taken to be 0.5. As usual, G and m p i stand for Newton's constant and Planck 
mass, respectively, and z(t) is the redshift at time t. Distances are measured in 
Mpc, where IMpc is 3.26 • 10 6 light years. t usually denotes the present time, and 
Tq stands for the present temperature. 
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2. Standard Cosmology: Successes and Problems 



The standard big bang cosmology rests on three theoretical pillars: the cosmo- 
logical principle, Einstein's general theory of relativity and a perfect fluid descrip- 
tion of matter. 

The cosmological principle 19 -* states that on large distance scales the Universe 
is homogeneous. From an observational point of view this is an extremely nontriv- 
ial statement. On small scales the Universe looks extremely inhomogeneous. The 
inhomogeneities of the solar system are obvious to everyone, and even by naked 
eye it is apparent that stars are not randomly distributed. They are bound into 
galaxies, dynamical entities whose visible radius is about 10 4 pc. Telescopic ob- 
servations show that galaxies are not randomly distributed, either. Dense clumps 
of galaxies can be identified as Abell clusters. In turn, Abell cluster positions are 
correlated to produce the large-scale structure dominated by sheets (or filaments) 
with typical scale 100 Mpc observed in recent redshift surveys 9 ). Until recently, 
every new survey probing the Universe to greater depth revealed new structures 
on the scale of the sample volume. In terms of the visible distribution of matter 
there was no evidence for large-scale homogeneity. This situation changed in 1992 
with the announcement 20 -* that a new redshift survey complete to a depth of about 
800 Mpc had discovered no prominent structures on scales larger than 100 Mpc. 
This is the first observational evidence from optical measurements in favor of the 
cosmological principle. However, to put this result in perspective we must keep in 
mind that the observed isotropy of the CMB temperature to better than 10~ 5 on 
large angular scales has been excellent evidence for the validity of the cosmological 
principle. 

The second theoretical pillar is general relativity, the theory which determines 
the dynamics of the Universe. According to general relativity, space-time is a 
smooth manifold. Together with the cosmological principle this tells us that it is 
possible to choose a family of hypersurfaces with maximal symmetry. These are 
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the homogeneous constant time hypersurfaces. The metric of these surfaces is 21 ) 



ds 2 = a(tf 



dr A 



+ r 2 {d6 2 + sin z 6 dip' z ) 



1 — kr 2 



(2.1) 



when using spherical polar coordinates. The constant A; is +1, or —1 for closed, 
flat or open surfaces respectively. The function a(t) is the scale factor of the 
Universe. By a coordinate choice, it could be set equal to 1 at any given time. 
However, the time dependence of a(t) indicates how the spatial sections evolve as 
a function of time. The full space-time metric is 



ds 2 = dt 2 -a(tf 



dr l 



+ r 2 (d6 2 +sin z 6 dtp' z ) 



1 — kr 2 



(2.2) 



According to the Einstein equivalence principle, test particles move on geodesies 
with respect to the metric given by ds 2 . This implies that the peculiar velocity v p 
obeys the equation 



Here, 



1 a F 



dx c 



(2.3) 



(2.4) 



x c being the comoving coordinates r, 9 and tp. Equation (2.3) implies that in an 
expanding Universe, the peculiar velocity v p decreases: 



vJt) ~ a-\t) 



(2.5) 



Therefore, trajectories with constant x c are geodesies and correspond to particles 
at rest. The velocity v p is the physical velocity relative to the expansion of the 
Universe (see Fig. 1). 
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Figure 1. Sketch of the expanding Universe. Concentric circles indicate space at 
fixed time, time increasing as the radius gets larger. Points at rest have constant 
comoving coordinates. Their world lines are straight lines through the origin (e.g. 
L). 

From observations 22 ^ it is known that the Universe is at present expanding. 
Looking at distant galaxies, we detect a redshift z of light which increases linearly 
with the distance d of the object: 



z(d) ~ d , 



(2.6) 



where the redshift z is defined by 



Ao 1 



(2.7) 



with Ao and A e being the wavelengths measured by us and by the source. The rela- 
tionship (2.6) is explained by taking galaxies to be at rest in comoving coordinates, 
and a(t) to be increasing. In this case for z <C 1 



H(t )d 



H(t) = -(t) 



(2.8) 



to being the present time. 

The most important consequence of general relativity for the history of the 
Universe is that it relates the expansion rate to the matter content. The Einstein 
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field equations follow from the action 



S = J d 4 xy/=g~(R- IQuGCm) (2.9) 

where R is the Ricci scalar curvature, g is the determinant of the metric, and Cm 
is the Lagrange density for matter. Evaluating the equations of motion obtained 
by varying (2.9) with respect to for a metric of the form (2.2) leads to the 
famous Friedmann-Robertson- Walker (FRW) equations 




(2.10) 



(, + %,). (2.11) 

a 6 

Note that for fixed energy density p, the evolution of a(t) depends in an important 
way on the pressure p. 

Besides explaining Hubble's redshift-distance relation (2.6), standard big bang 
cosmology makes two key quantitative predictions: the existence of a black body 
cosmic microwave background 23 -* and nucleosynthesis, the generation of light elements 23-25 ) .| 

Consider ordinary matter made up of atoms in an expanding Universe. The en- 
ergy density in matter scales as a(t)~ 3 , and the temperature T(t) as a _1 (t). Thus, 
standard cosmology predicts that as we go back in time, the Universe was warmer. 
In particular, at a critical temperature T rec , matter becomes ionized. Before t rec 
(the time corresponding to T rec ) the Universe was opaque to photons, after t rec it 
was transparent. To be more precise, T rec is the temperature when photons fall out 
of thermal equilibrium. Thereafter they propagate without scattering. The black 
body nature of the spectrum of photons is maintained, but the temperature red- 
shifts. Hence, the big bang model predicts 23 ' 26 ) a black body spectrum of photons 
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of temperature 

T — T rec z~g C (2-12) 
where the cosmological redshift z(t) is given by 

1 + *(*) = ^ (2.13) 
a(t) 

and z rec = z(t rec ) ~ 10 3 . 

In 1965, Penzias and Wilson 27 ) discovered this remnant black body radiation 
at a temperature of about 3° K. Since the spectrum peaks in the microwave region 
it is now called CMB (cosmic microwave background). Recent satellite (COBE) 28 ) 
and rocket 29 ) experiments have confirmed the black body nature of the CMB to 
very high accuracy. The temperature is 2.73° K = Tq. 

Given the existence of the CMB, we know that matter has two components: 
dust (with energy density p m (t)) and radiation (with density p r (t)). At the present 
time to, p m (t) ^> p r (t). The radiation energy density is determined by To, and the 
matter energy density can be estimated by analyzing the dynamics of galaxies and 
clusters and using the virial theorem. However, since 

p m (t) ~ a(t)~ 3 

(2.14) 

p r (t) ~ a(t) 4 , 

as we go back in time the fraction of energy density in radiation increases, and the 
two components become equal at t eq , the time of equal matter and radiation. The 
corresponding redshift is 

z eq = nh 5 *W 4 (2.15) 

where 

^ = — (t ), (2.16) 

Per 

Per being the density for a spatially flat Universe (the critical density), and h§o in 
the value of H in units of 50 km s _1 Mpc -1 . 
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The time t eq is important for structure formation 30 ). It is only after t eq that 
perturbations on scales smaller than the Hubble radius H^it) can grow. Before 
then, the radiation pressure prevents growth. A temperature-time plot of the early 
Universe is sketched in Fig. 2. Note that t eq < tree- 



Figure 2. Temperature-time diagram of standard big bang cosmology. The present 
time, time of decoupling and time of equal matter and radiation are to, t rec and 
t eq , respectively. The Universe is radiation dominated before t eq (Region A), and 
matter dominated in Region B. Before and after t rec the Universe is opaque and 
transparent, respectively, to microwave photons. 



The second quantitative prediction of standard big bang cosmology concerns 
nucleosynthesis. Above a temperature of about 10 9 ° K, the nuclear interactions are 
sufficiently fast to prevent neutrons and protons from fusing. However, below that 
temperature, it is thermodynamically favorable 23-25 -* for neutrons and protons to 
fuse and form deuterium, helium 3, helium 4 and lithium 7 through a long and 
interconnected chain of reactions. The resulting light element abundances depend 
sensitively on the expansion rate of the Universe and onOg, the fraction of energy 
density ps at present in baryons relative to the critical density p c . In Fig. 3, 
recent 31 ^ theoretical calculations of the abundances are shown and compared with 
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observations. Demanding agreement with all abundances leaves only a narrow 
window 



Figure 3. Light element abundances as a function of the baryon to entropy ratio 
rj. Y p is the 4 He mass fraction, calculated assuming three light neutrino species, 
for two different values of the neutrino half life (in minutes). The abundances 
of £>+ 3 He and of 7 Li are shown as ratios of their number density relative to the 
number density of H. The horizontal lines indicate limits from observations: upper 
limit of 10 -4 for the .D+ 3 He abundance, and upper limits on 7 Li from observations 
of dwarf stars. The 7 Li curve is shown with ±2<r errors. Combining 4 He and 3 He 
+D limits leaves only a small window for rj which is allowed. A major success of 
primordial nucleosynthesis is that the 7 Li abundance matches well. Note that t]iq is 
rj in units of 1CT 10 . 
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3 x 1(T 1U < rj < 10 



,-9 



(2.17) 



where 77 is the ratio of baryon number density to entropy density s 




(2.18) 



From (2.17), it follows that Qb is constrained: 



0.01 < n B h 2 < 0.035. 



(2.19) 



In particular, if the Universe is spatially flat, there must be nonbaryonic dark 
matter. We will return to the dark matter issue shortly. 

In summary, the three observational pillars of standard big bang cosmology are 
Hubble's redshift- distance relation, the existence and black body spectrum of the 
CMB, and the concordance between observed and theoretically determined light 
element abundances. 

Standard big bang cosmology is faced with several important problems: the 
age, dark matter, homogeneity, flatness, and formation of structure problems (this 
list is not exhaustive 32 ) , but contains what I consider to be the key problems). 
In addition, standard big bang cosmology does not explain the small value of the 
cosmological constant - a problem which no present cosmological model addresses 
(see, however, Ref. 33 for some recent progress on this issue). Of the above five 
problems, only the first is a possible conflict of the theory with observations. The 
other four are questions which are left unanswered by the theory. Extensions of 
the standard model are needed to address these issues. 

Globular cluster ages have been estimated to lie in the range 13 — 18 x 10 9 
years 34 ). Nuclear cosmochronology gives an age of structures in the range (10 — 
20) x 10 9 yr. However, big bang cosmology (in the absence of a cosmological 



12 



constant) predicts an age 

7 

r=--10 2 yr- (2.20) 

Thus, theory and observations are only compatible if h < 0.55. 

The dark matter problem 35 ) problem has various aspects. There is more mat- 
ter in galaxies than is visible in stars. This follows by studying galactic velocity 
rotation curves and observing that the velocity remains constant beyond the visible 
radius of the galaxy. Whereas 34 ) the contribution of stars to Q is less than 0.01, 
galaxies as a whole contribute > 0.02 to VI. The second level of the dark matter 
problem lies in clusters. By studying cluster dynamics it can be inferred that the 
contribution to Q exceeds 0.1. Finally, observations on the largest scales give an 
even larger contribution to Q. From Virgo infall it follows that Q > 0.3, and from 
large scale velocity measurements (POTENT) 36 ) or from infrared galaxy (IRAS) 
surveys 37 ) it follows that 

0.5<n<3. (2.21) 

Inflation predicts Q ~ 1, unless fine tuned initial conditions are chosen 38 ). 

The bottom line of the dark matter problem is that more mass is observed 
than can be in light (hence the name "dark matter"), and that most of this missing 
matter must be nonbaryonic (since by (2.19) Qb < 0.14). Note that there must 
be some baryonic dark matter, and that all of the galactic dark matter could 
be baryonic. This comment will be relevant when discussing cosmic strings and 
galaxy formation in Section 6. If h = 0.5, then most of the cluster dark matter 
must be nonbaryonic. Standard cosmology does not address the issue of what the 
dark matter is. New cosmology can provide a solution. Quantum field theory 
models of matter in the Universe give rise to several candidate particles which 
could constitute the dark matter 39 ). 

The last three problems are the classic problems which the inflationary Uni- 
verse scenario addresses 2 ). In Fig. 4, the homogeneity (or horizon) problem is 
illustrated. As is sketched, the region £p(t rec ) over which the CMB is observed to 
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be homogeneous to better than one part in 10 5 is much larger than the forward 
light cone £f(t rec ) at t rec , which is the maximal distance over which microphysical 
forces could have caused the homogeneity: 

£ p (trec) = a(t rec ) j dt <r\t) ~ 3t ^1 - (^-^j j a(t rec ) (2.22) 



if (tree) = a(t r ec) j dt a (t) ~ 3t Q ' t r ' ec a(t rec ) . (2.23) 


From the above equations it is obvious that £ p (t eq ) 3> £f(t eq ). Hence, standard 
cosmology cannot explain the observed isotropy of the CMB. 



Figure 4. A sketch of the homogeneity problem: the past light cone £ p (t) at the 
time t rec of last scattering is much larger than the forward light cone £f(t) at t rec . 



In standard cosmology and in an expanding Universe, Q = 1 is an unstable 
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fixed point. This can be seen as follows. For a spatially flat Universe (Q = 1) 

ff 2 = ^Pc, (2.24) 

whereas for a nonflat Universe 

# 2 + eT 2 =^p, (2.25) 

with 

E =(^r < 2 - 26 > 

The quantity e is proportional to s~ 2//3 , where s is the entropy density. Hence, in 
standard cosmology, e is constant. Combining (2.24) and (2.25) gives 

P^ = ^ £ -H~T-\ (2.27) 
Pc SnG Pc V ; 

Thus, as the temperature decreases, Q — l increases. In fact, in order to explain the 
present small value of fl — 1 ~ 0(1), the initial energy density had to be extremely 
close to critical density. For example, at T = 10 15 GeV, (2.27) implies 

P -^- ~ 1(T 50 . (2.28) 

Pc 

What is the origin of these fine tuned initial conditions? This is the flatness problem 
of standard cosmology. 

The last problem of the standard cosmological model I will mention is the "for- 
mation of structure problem." Observations indicate that galaxies and even clusters 
of galaxies have nonrandom correlations on scales larger than 50 Mpc 9 ' 40 ' 41 ). This 
scale is comparable to the comoving horizon at t eq . Thus, if the initial density per- 
turbations were produced much before t eq , the correlations cannot be explained by 
a causal mechanism. Gravity alone is, in general, too weak to build up correlations 
on the scale of clusters after t eq (see, however, the explosion scenario of Ref. 42). 
Hence, the two questions of what generates the primordial density perturbations 
and what causes the observed correlations, do not have an answer in the context 
of standard cosmology. This problem is illustrated by Fig. 5. 
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Figure 5. A sketch (conformal separation vs. time) of the formation of structure 
problem: the comoving separation d c between two clusters is larger than the for- 
ward light cone at time t rec - 

In 1981, based on previous work by many other people (see e.g., Refs. 43-45 
for a detailed bibliography of early work), Guth realized 2 ) that having a suffi- 
ciently long phase in the very early Universe during which the scale factor expands 
exponentially 

a(t) ~ e m (2.29) 

can potentially solve the three last problems listed above. This phase of exponential 
expansion is called the de Sitter or inflationary phase. 

In Fig. 6 it is sketched how a period of inflation can solve the homogeneity 
problem. U shall denote the onset of inflation, tR the end. At — tR — U is the 
period of inflation. During inflation, the forward light cone increases exponentially 
compared to a model without inflation, whereas the past light cone is not affected 
for t > tR. Hence, provided At is sufficiently large, £f(t rec ) will be greater than 
£p(trec)- The condition on At depends on the temperature Tr corresponding to 
time tR, the temperature of reheating. Demanding that if(tR) > £ p (tR) we find, 
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Figure 6. A sketch (physical coordinates vs. time) of the solution of the homo- 
geneity problem. Due to exponential expansion the forward light cone is larger 
than the past light cone at t Tec . The dashed line is the apparent horizon or Hubble 
radius. 

using (2.22) and (2.23), the following criterion 

1/2 



e AtH > 



W^^i'A =f^)~ 10 27 (2.30) 
if(t R ) \t R J \T J 



for Tji ~ 10 14 GeV and To ~ 10~ 13 GeV (the present microwave background tem- 
perature). Thus, in order to solve the homogeneity problem, a period of inflation 
with 

At>50H- 1 (2.31) 

is required. 
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Inflation also can solve the flatness problem 2 ' 46 ) . The key point is that the 
entropy density s is no longer constant. As will be explained later, the temperatures 
at U and tji are essentially equal. Hence, the entropy increases during inflation 
by a factor exp(3iJAi). Thus, e decreases by a factor of exp(—2HAt). With the 
numbers used in (2.30): 

Cafter ~ 10~ 54 e be fore ■ (2.32) 

Hence, (p — p c )/p can be of order 1 both at U and at the present time. In fact, if 
inflation occurs at all, the theory then predicts that at the present time Q = 1 to 
a high accuracy (now Q < 1 would require special initial conditions). 

What was said above can be rephrased geometrically: during inflation, the 
curvature radius of the Universe - measured on a fixed physical scale - increases 
exponentially. Thus, a piece of space looks essentially flat after inflation even if it 
had measurable curvature before. 

Most importantly, inflation provides a mechanism which in a causal way gen- 
erates the primordial perturbations required for galaxies, clusters and even larger 
objects 47-50 ). In inflationary Universe models, the "apparent" horizon 3t and the 
"actual" horizon (the forward light cone) do not coincide at late times. Provided 
(2.30) is satisfied, then (as sketched in Fig. 7) all scales within our apparent hori- 
zon were inside the actual horizon since U. Thus, it is in principle possible to have 
a causal generation mechanism for perturbations. 

As will be shown at the beginning of the following chapter, in order to obtain 
inflation it is necessary to describe matter in terms of field theory. In particular, 
there must be a phase transition in the early Universe during which the equation of 
state of matter changes, leading to different expansion rates of the Universe before 
and after the transition. 

In order to obtain sufficient inflation 2 \ the phase transition must be slow. 
For particle physics motivated matter theories this is hard to achieve as will again 
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be shown in the next chapter. However, matter theories generically predict phase 
transitions, and the latter lead to the formation of topological defects. 

Topological defect models provide a second mechanism for producing structures 
in the Universe on large scales. To be specific, consider a model giving rise to 
cosmic strings (see Chapter 4). The strings produced during the phase transition 
will form a network of random walks. This network has infinite extent and gives 
rise to nonrandom correlations on all scales for point objects seeded by the strings 
(see Chapter 6). This and similar topological defect models for structure formation 
will be the focus of a substantial part of these lectures. 



Figure 7. A sketch (physical coordinates vs. time) of the solution of the forma- 
tion of structure problem. The separation d c between two clusters is always smaller 
than the forward light cone. The dashed line is the Hubble radius iJ _1 (t). 
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3. Phase Transitions and Topological Defects 



3.1. Scalar Fields and Cosmology 

Let us for a moment return to the inflationary Universe. As stated in (2.29), 
the requirement for inflation is to have a time interval during which the scale factor 
a(t) expands exponentially. Recall form (2.10) and (2.11) that for a spatially flat 
FRW Universe the Einstein equations reduce to 



a\ 2 8ttG 



a J 3 



P (3-1) 



and 

2 



2 - + ( _ ] = SnGp (3.2) 
a \a 



However, using a(t) = e tH , the left hand side of (3.2) can be evaluated using (3.1) 
to give 



2 

a fa 



2- + - = 8nGp. (3.3) 
a \a J 

In order that (3.2) and (3.3) agree, the equation of state of matter must be 

P=-P (3.4) 

Therefore, 'new cosmology' (i.e. description of matter in terms of fields) is required 
in order to obtain inflation. 

It is not hard to show that an equation of state like (3.4) with negative pressure 
can be obtained if matter is described in terms of scalar fields. Consider the 
Lagrangian L(ip) for a theory of a scalar field f(x,t): 

C(<p) = ±d ll <pd' t <p-V(<p). (3.5) 
Given the Lagrangian, the energy-momentum tensor can be determined as in 
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any classical field theory (see eq. Ref. 51). In a Universe with FRW metric 

= diag (l, -a 2 (t), -a 2 (t), -a 2 {t)) (3.6) 
we obtain for p = Tqq and p — ^ Tu: 

p(x,t) = l -(p 2 (x,i) + V 2 (VV^) 2 + W 

3 (3-7) 
p(x,t) = ^X>* = ^ 2 (x,t)- V 2 (vy?) 2 -U(y?). 
i=l 

Thus, if </?(x, tj) =const and 0(x, ti) = at some initial time U and V(ip(x, ti)) > 0, 
then the equation of state becomes p = —p and leads to inflation. 



Figure 8. A sketch of two potentials which can give rise to inflation. 

Two examples which give inflation are shown in Fig. 8. In (a), inflation occurs 
at the stable fixed point (p(x,t{) = = (p(x,ti). However, this model is ruled out 
by observation: the inflationary phase has no ending. V(0) acts as a permanent 
nonvanishing cosmological constant. In (b), a finite period of inflation can arise 
if v?(x) is trapped at the local minimum tp — with y?(x) = 0. However, in this 
case </?(x) can make a sudden transition at some time tR > t{ through the potential 
barrier and move to y?(x) = a. Thus, for ti < t < tR the Universe expands 
exponentially, whereas for t > tR the contribution of if to the expansion of the 
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Universe vanishes and we get the usual FRW cosmology. There are two obvious 
questions: how does the transition occur and why should the scalar field have 
V(tp) = at the global minimum. In the following section the first question will 
be addressed. The second question is part of the cosmological constant problem 
for which there is as yet no convincing explanation. Before studying the dynamics 
of the phase transition, we need to digress and discuss finite temperature effects. 

3.2. Finite Temperature Field Theory 

The evolution of particles in vacuum and in a thermal bath are very different. 
Similarly, the evolution of fields changes when coupled to a thermal bath. Under 
certain conditions, the changes may be absorbed in a temperature dependent po- 
tential, the finite temperature effective potential 52 ) . Here, a heuristic derivation 
of this potential will be given. The reader is referred to Ref.53 or to the original 
articles 52 ^ for the actual derivation. 

We assume that the scalar field <p(x,t) is coupled to a thermal bath which is 
represented by a second scalar field ip(x, t) in thermal equilibrium. The Lagrangian 
for ip is 

C = \d^d^ - V(tp) - , (3.8) 

where A is a coupling constant. The action from which the equations of motion 
are derived is 

S = J d 4 x (3.9) 

where g is the determinant of the metric (3.6). The resulting equation of motion 
for </?(x, t) is 

y3 + 3#yj -a" 2 sj 2 ip = -V'(<p) - A> V (3.10) 
If ip is in thermal equilibrium, we may replace ^ 2 by its thermal expectation value 
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< ip 2 > T . Now, 



< ip 2 > r ~ T 2 



(3-11) 



which can be seen as follows: in thermal equilibrium, the energy density of ip equals 
that of one degree of freedom in the thermal bath. In particular, the potential 
energy density V(ip) of ip is of that order of magnitude. Let 

V(iP) = A^ 4 (3.12) 

with a coupling constant A^, which we take to be of the order 1 (if \^ is too 
small, ip will not be in thermal equilibrium). Since the thermal energy density is 
proportional to T 4 , (3.11) follows. (3.10) can be rewritten as 

ip + 3H(p> - a" 2 v 2 V? = -V>r((p), (3.13) 

where 

V T (<p) = V(<p) + l -\T 2 V 2 (3.14) 

is called the finite temperature effective potential. Note that in (3.14), A has been 
rescaled to absorb the constant of proportionality in (3.11). 

These considerations will now be applied to Example A, a scalar field model 
with potential 

V{<p) = \\(<P 2 -V 2 ) 2 (3.15) 

(77 is called the scale of symmetry breaking). The finite temperature effective 
potential becomes (see Fig. 9) 

Vritp) = ~ \ (W ~ AT 2 ) V 2 + \ \ V 4 . (3.16) 

For very high temperatures, the effective mass term is positive and hence the 
energetically favorable state is < (p >= 0. For very low temperatures, on the other 
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hand, the mass term has a negative sign which leads to spontaneous symmetry 
breaking. The temperature at which the mass term vanishes defines the critical 
temperature T c 

T c = \~ 1 / 2 \ 1 / 2 i] . (3.17) 



Figure 9. The finite temperature effective potential for Example A. 

As Example B, consider a theory with potential 

V(<p) = \^ - \ (a + % 3 + l -ab V 2 (3.18) 

with > b > 0. The finite temperature effective potential is obtained by adding 
\\T 2 ip 2 to the right hand side of (3.18). Vr((p) is sketched in Fig. 10 for various 
values of T. The critical temperature T c is defined as the temperature when the 
two minima of Vr((p) become degenerate. 

It is important to note that the use of finite temperature effective potential 
methods is only legitimate if the system is in thermal equilibrium. This point was 
stressed in Refs. 54 and 55, although the fact should be obvious from the derivation 
given above. To be more precise, we require the ip field to be in thermal equilib- 
rium and the coupling constant A of (3.8) which mediates the energy exchange 
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Figure 10. The finite temperature effective potential for Example B. 

between the if and ip fields to be large. However, (see e.g. Ref. 18), observational 
constraints stemming from the amplitude of the primordial energy density fluctua- 
tion spectrum force the self coupling constant A of if to be extremely small. Since 
at one loop order, the interaction term ^Xip 2 ^ 2 induces contributions to A, it is 
unnatural to have A very small and A unsuppressed. Hence, in many inflationary 
Universe models - in particular in new inflation 56 ) and in chaotic inflation 54 ) - finite 
temperature effective potential methods are inapplicable. 

3.3. Phase Transitions 

The temperature dependence of the finite temperature effective potential in 
quantum field theory leads to phase transitions in the very early Universe. These 
transitions are either first or second order. 

Example A of the previous section provides a model in which the transition is 
second order (see Fig. 9). For T ^> T c , the expectation value of the scalar field if 
vanishes at all points x in space: 

<y?(x)>=0. (3.19) 

For T < T c , this value of < ip(x) > becomes unstable and < </?(x) > evolves 
smoothly in time to a new value ±i]. The direction is determined by thermal 
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and quantum fluctuations and is therefore not uniform in space. There will be 
domains of average radius £(£) in which < </?(x) > is coherent. By causality, the 
coherence length is bounded from above by the horizon. However, typical values of 
are proportional to A -1 ?? -1 if <p was in thermal equilibrium before the phase 
transition 1 ^ . 

In condensed matter physics, a transition of the above type is said to proceed 
by spinodal decomposition 57 ) , triggered by a rapid quench. 

In Example B of the previous section, (see Fig. 10) the phase transition is first 
order. For T > T c , the expectation value < (p(x) > is approximately 0, the min- 
imum of the high temperature effective potential. Provided the zero temperature 
potential has a sufficiently high barrier separating the metastable state ip — from 
the global minimum (compared to the energy density in thermal fluctuations at 
T = T c ), then </?(x) will remain trapped at (p — also for T <T C . In the notation 
of Ref. 58, the field ip is trapped in the false vacuum. After some time (determined 
again by the potential barrier), the false vacuum will decay by quantum tunnelling. 

Tunnelling in quantum field theory was discussed in Refs. 58-61 (for reviews 
see e.g., Refs. 62 and 53). The transition proceeds by bubble nucleation. There is 
a probability per unit time and volume that at a point x in space a bubble of "true 
vacuum" </?(x) = a will nucleate. The nucleation radius is microscopical. As long 
as the potential barrier is large, the bubble radius will increase with the speed of 
light after nucleation. Thus, a bubble of <p — a expands in a surrounding "sea" of 
false vacuum ip — 0. 

To conclude, let us stress the most important differences between the two types 
of phase transitions discussed above. In a second order transition, the dynamics 
is determined mainly by classical physics. The transition occurs homogeneously 
in space (apart from the phase boundaries which - as discussed below - become 
topological defects), and < <p(x) > evolves continuously in time. In first order tran- 
sitions, quantum mechanics is essential. The process is extremely inhomogeneous, 
and < (p(x) > is discontinuous as a function of time. 
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3.4. Particle physics connection 

Up to this point in these lectures, I have only discussed scalar field toy models. 
The actual case of interest is a unified gauge theory of strong, weak and electro- 
magnetic interactions. Such a theory (see e.g. Ref. 63) contains massless fermion 
fields ip, gauge fields and scalar fields tp. The Lagrangian L(ip, A^, tp) is in- 
variant under the action of some internal symmetry group G. This grand unified 
symmetry is broken spontaneously (see e.g. Ref. 64) in one or several stages for 
the 'standard model' symmetry group SU(3) x SU(2) x U(l): 

G -> H -> . . . -> SU(3) x SU(2) x 17(1). (3.20) 

Spontaneous symmetry breaking is achieved by postulating that a scalar (Higgs) 
field takes on an expectation value which is not invariant under the full group G, but 
only under the subgroup H. In order that this be possible, the potential V(ip) for 
this Higgs field <p must have a global minimum at <p ^ but a high temperature 
minimum at <p — 0. Therefore, such a model will lead to a symmetry breaking 
phase transition in the early Universe. Hence, the possibility of topological defect 
formation arises. 

As a toy model, consider a theory of a two component real scalar field 




(3.21) 



with Lagrangian 

C{p)= l -d^p-V{^) ) (3.22) 
and symmetry breaking potential 

V( V )= 1 -\( V 2 - V 2 ) 2 . (3.23) 

Summation over the internal index % — 1, 2 is implicit in (3.22) and (3.23). Obvi- 
ously, the Lagrangian is invariant under a U(l) rotation U(a) in an internal field 
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space 3? . However, a choice of vacuum state 



Vvac = V 



( 







1 



) 



(3.24) 



breaks this symmetry completely. 

The concept of vacuum manifold Ai will be crucial in the following sections. 
In general, M. is defined is the set of field configurations which minimize the 
free energy modular gauge transformations. For a gauge theory with symmetry 
breaking scheme as in the first stage of (3.20) 



I will now argue that in 'realistic' particle physics models, a symmetry breaking 
phase transition proceeds too fast to allow for inflation. Hence, the topological 
defects which form in the phase transition are not inflated away (i.e. moved beyond 
the present Hubble radius by the exponential expansion of the forward 

light cone), but remain inside the Hubble radius and are therefore important for 
cosmology. 

The argument is based on the equation of motion for the Higgs field. Neglect- 
ing inhomogeneitis in if, interactions with other fields, and the expansion of the 
Universe, this equation becomes 



M ~ G/H. 



(3.25) 



In the above toy model, 



M = W : V(<f) = V min = 0}, 



(3.26) 



if ~ — — — ~ Xf] if 

dip 



(3.27) 



The typical time scale r for the solution is 



r = A-V. 



(3.28) 
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Sufficient inflation requires 

r^H- 1 . (3.29) 
However, from the FRW equation (3.1) it follows that 

H 2 - ^V(O) = |a<V, (3.30) 

and hence 

r^H- 1 (3.31) 

unless 7] > rripi. 

In order to obtain inflation, it is in general 54 > 65 ) necessary to introduce a scalar 
field (f which is not the Higgs field of a unified field theory. We shall not further 
consider such models and restrict our attention to matter theories motivated by 
particle physics 

3.5. Classification of Topological Defects 



Different particle physics models admit different types of topological defects. 
The topology of the vacuum manifold M. determines the type of defect. 

The classification of defects is based on homotopy classes of M. To introduce 
the ideas (see e.g. Ref. 66), consider maps ip from S n to M, n being an integer. 
Two maps ip\ and ip2 are called homotopically equivalent (ipi ~ ^2) if there exists 
a continuous one parameter family of maps 

i){t) :S n ^M t e [0, 1] (3.32) 

with 

V>(0) = Vi , ^(1) = V>2- (3.33) 

The n'th homotopy 'group' of M, U n (Ai), is the set of all homotopy classes of 
maps S n — > M.. Except for n = 0, Tl n (M) is a group. 
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Let us consider a symmetry breaking pattern 

G^H (3.34) 

of a group G with 

U (G) = ni(G) = 1 (3.35) 

to a subgroup ff. As stated in (3.25), the vacuum manifold M is isomorphic to 
G/H. Using exact sequences 66 ) it can be shown X ' 4 ) that 

IIi(.M)=n (#) (3.36) 

and 

II 2 (A<) =Il 1 (H), (3.37) 

assuming 112(G) = 1 for the second relation to hold. 

Let us consider two relevant examples. The smallest grand unified group is 
£77(5). If we consider the symmetry breaking pattern 

£77(5) -> SU(3) x SU(2) x U(l), (3.38) 

we can use (3.36) and (3.37) to conclude that 

Ui(M) = 1 (3.39) 

and 

IL 2 (M) = IL 1 {U(1))=Z. (3.40) 

As will be shown in the following section, (3.39) and (3.40) imply that no cosmic 
strings but monopoles form during the phase transition associated with (3.38). 
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As a second example, consider G = 5*0(10) and the symmetry breaking pattern 

50(10) -> 517(5) x Z 2 . (3.41) 

In this example 

U l (M) = Z 2 , (3.42) 

and hence cosmic strings form during the symmetry breaking phase transition. 

Let us give a brief overview of the classification of topological defects. To be 
specific, consider a theory with an n-component real scalar field ip with symmetry 
breaking potential (3.23). 

There are various types of local and global topological defects 3 ) (regions of 
trapped energy density) depending on the number of components of ip. The words 
"local" and "global" refer to whether the symmetry which is broken is a gauge or 
global symmetry. In the case of local symmetries, the topological defects have a well 
defined core outside of which p> contains no energy density in spite of nonvanishing 
gradients Vy?: the gauge fields can absorb the gradient, i.e., D^ip = when 
djjip 7^ 0, where the covariant derivative is defined by 

D fl = d fl -igAp, (3.43) 

g being the gauge coupling constant. Global topological defects, however, have 
long range density fields and forces. 

Table 1 contains a list of topological defects with their topological character- 
istic. A "v" marks acceptable theories, a "x" theories which are in conflict with 
observations (for rj ~ 10 16 GeV). 
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In the following sections of this chapter we will construct the various types of 
defects. The cosmological implications will be discussed in later chapters. 

The easiest defect to construct is the domain wall. Consider n — 1 (or more 
generally a theory with IIq(M) ^ 1). In this case, the vacuum manifold consists 
of two points 

ip vac = ±V- (3.44) 

During the symmetry breaking phase transition, regions in physical space 1Z 3 with 
ip = ±rj will form. These regions are separated by two dimensional surfaces (walls) 
with tp M. (see Fig. 11). These are the domain walls. Since p ^ M. in the 
walls, V(tp) > and hence the walls carry energy per unit area. Via the usual 
gravitational force, this energy can act as a seed for structures in the Universe. 



Figure 11. A two dimensional cross section through space showing a domain wall 
(DW) separating a region with ip — r] (+) from a neighboring region with ip = —r\ 

(")• 



3.6. Cosmic Strings 

Consider a theory in which matter consists of a gauge field and a complex 
scalar field <fi whose dynamics is given by the Lagrangian 

C= l - D^D^ - V{(f>) + \ F^v (3.45) 

where F^ v is the field strength tensor. The potential V(<p) has the symmetry 
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breaking "Mexican hat" shape (see Figure 12): 



v(<t>) = l -\(\<t ) \ 2 - v 2 ) 2 . 



(3.46) 



Hence, the vacuum manifold Ai, the space of minimum energy density configura- 
tions, is a circle S 1 . 



Figure 12. The zero temperature potential energy of the complex scalar field used 
in the cosmic string model. 

The theory described by (3.45) and (3.46) admits one dimensional topological 
defects, cosmic strings. In the Abelian Higgs model of this example the string 
solutions were first found by Nielsen and Olesen 67 -*. It is possible to construct 
string configurations which are translationally invariant along the z axis. On a 
circle C in the x — y plane with radius r(see Fig. 13), the boundary conditions for 
are 



0(r, 9) = r] e 



(3.47) 



where 6 is the polar angle along C. 
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Figure 13. Sketch of the cosmic string construction of section 3.6 . See text for 
notation. 

The configuration (3.47) has winding number 1: at all points of the circle, <fi 
takes on values in Ai, and as tp varies from to 2it, <fi winds once round Ai. By 
continuity it follows that there must be a point p on the disk D bounded by C where 
= 0. By translational symmetry there is a line of points with = 0. This line 
is the center of the cosmic string. The cosmic string is a line of trapped potential 
energy. In order to minimize the total energy given the prescribed topology (i.e., 
winding number), the thickness of the string (i.e., radius over which V(<f>) deviates 
significantly from 0) must be finite. As first shown in Ref. 67, the width w of a 
string is 

w ~ A^V 1 , (3-48) 
from which it follows that the mass per unit length \i is 

fi ~ rf , (3.49) 

i.e., independent of the self-coupling constant A. 

Cosmic strings arise in any model in which the vacuum manifold satisfies the 
topological criterion 

Iii(M)^l. (3.50) 
Any field configuration 0(x) is characterized by an integer n, the element of Ui(Ai) 
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corresponding to 0(x). (Roughly speaking, n can be viewed as the number of times 
the map ip from C to M. covers M. ). 

A cosmic string is an example of a topological defect . A topological defect has 
a well-defined core, a region in space where ^ M. and hence V{fa) > 0. There 
is an associated winding number, and it is quantized. Hence, a topological defect 
is stable. Furthermore, topological defects exist for theories with global and local 
symmetry groups. 

3.7. Monopoles 

If the theory contains three real scalar fields fa with potential (3.46) (if |0| 2 = 

3 

<fij), then n.2(.M) 7^ 1 and monopoles result. The construction of a monopole 

i=l 

configuration is illustrated in Fig. 14. As the origin in physical space we select 
a point which is to become the center of the monopole. Consider a sphere S r of 
radius r surrounding this point. A spherically symmetric monopole configuration 
is obtained by the identity map 

S r -> M = S 2 

(r,e,<p)^>(0,<p). (3 ' 51) 

This configuration has winding number 1. Since the winding number of maps 
S 2 — > S 2 is quantized, it cannot change as r varies. Thus, the only way to obtain 
a single valued field configuration at r = is for tp(r, 8, tp) to leave M. as r — > 0. 
In particular, there is a point (e.g., r = 0) for which = 0. This is the center of 
the monopole. We see that monopoles are topological defects: they contain a core, 
have quantized winding number and are stable. 
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Figure 14. Construction of a monopole: left is physical space, right the vacuum 
manifold. The field configuration maps spheres in space onto M.. However, a 
core region of space near the origin is mapped onto field values not in M.. 

3.8. Global Textures 

Next, consider a theory of four real scalar fields given by the Lagrangian 

C = \d^cf> - V(cf>) (3.52) 

with 

In this case, the vacuum manifold is M. = S s with topology 

U 3 (M)^1, (3.54) 
and the corresponding defects are the global textures 1 ' 68 ' 69 ). 
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Figure 15. Construction of a global texture: left is physical space, right the vacuum 
manifold. The field configuration is a map from space to the vacuum manifold (see 
text). 

Textures, however, are quite different than the previous topological defects. 
The texture construction will render this manifest (Fig. 15). To construct a 
radially symmetric texture, we give a field configuration 4>{x) which maps physical 
space onto M.. The origin in space (an arbitrary point which will be the center 
of the texture) is mapped onto the north pole N of M. Spheres surrounding 
are mapped onto spheres surrounding N. In particular, some sphere with radius 
r c (t) is mapped onto the equator sphere of M. The distance r c (t) can be defined 
as the radius of the texture. Inside this sphere, <f>(x) covers half the vacuum 
manifold. Finally, the sphere at infinity is mapped onto the south pole of M.. The 
configuration 0(x) can be parameterized by 69 ) 

(f)(x,y,z) = cosx(r), sinx(r)-, sinx(r)-, sinx(r)-) (3.55) 
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in terms of a function x( r ) with x(0) = an d x(°°) = 71 ■ Note that at all points 
in space, 0(x) lies in Ai. There is no defect core. All the energy is spatial gradient 
(and possibly kinetic) energy. 

In a cosmological context, there is infinite energy available in an infinite space. 
Hence, it is not necessary that x( r ) ^ vr as r ^ oo. We can have 

Xir) Xmax < 7T as T -> oo . (3.56) 
In this case, only a fraction 

n = ^ - Sin2Xmax (3.57) 

7T 27T V 7 

of the vacuum manifold is covered: the winding number n is not quantized. This 
is a reflection of the fact that whereas topologically nontrivial maps from S* 3 to S 3 
exist, all maps from i? 3 to S* 3 can be deformed to the trivial map. 

Textures in i? 3 are unstable. For the configuration described above, the insta- 
bility means that r c (t) — > as t increases: the texture collapses. When r c (t) is 
microscopical, there will be sufficient energy inside the core to cause 0(0) to leave 
M., pass through and equilibrate at x(0) = vr: the texture unwinds. 

A further difference compared to topological defects: textures are relevant only 
for theories with global symmetry. Since all the energy is in spatial gradients, for a 
local theory the gauge fields can reorient themselves such as to cancel the energy: 

D^ = 0. (3.58) 

Therefore, it is reasonable to regard textures as an example of a new class 
of defects, semitopological defects. In contrast to topological defects, there is no 
core, and (f>(x)eA4 for all x- In particular, there is no potential energy. Second, 
the winding number is not quantized, and hence the defects are unstable. Finally, 
they exist only in theories with a global internal symmetry. 
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4. Formation and Evolution of Topological Defects 



4.1. Kibble Mechanism 

The Kibble mechanism 1 -' ensures that in theories which admit topological or 
semitopological defects, such defects will be produced during a phase transition in 
the very early Universe. 

Consider a mechanical toy model, first introduced by Mazenko, Unruh and 
Wald 55 ) in the context of inflationary Universe models, which is useful in under- 
standing the scalar field evolution. Consider (see Fig. 16) a lattice of points on a 
flat table. At each point, a pencil is pivoted. It is free to rotate and oscillate. The 
tips of nearest neighbor pencils are connected with springs (to mimic the spatial 
gradient terms in the scalar field Lagrangian). Newtonian gravity creates a poten- 
tial energy V((p) for each pencil (</? is the angle relative to the vertical direction). 
V(ip) is minimized for \ip\ = r\ (in our toy model r\ = ir/2). Hence, the Lagrangian 
of this pencil model is analogous to that of a scalar field with symmetry breaking 
potential (3.46). 



Figure 16. The pencil model: the potential energy of a simple pencil has the same 
form as that of scalar fields used for spontaneous symmetry breaking. The springs 
connecting nearest neighbor pencils give rise to contributions to the energy which 
mimic spatial gradient terms in field theory. 
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At high temperatures T ^> T c , all pencils undergo large amplitude high fre- 
quency oscillations. However, by causality, the phases of oscillation of pencils with 
large separation s are uncorrelated. For a system in thermal equilibrium, the length 
s beyond which phases are random is the correlation length £(£). By causality there 
is an apriori causality bound on £: 



where t c is the causal horizon, at temperature T c . 

The critical temperature T c is the temperature at which the thermal energy 
is equal to the energy a pencil needs to jump from horizontal to vertical position. 
For T < T c , all pencils want to lie flat on the table. However, their orientations 
are random beyond a distance of £(t c ). 

The boundaries between the domains of correlated orientation become topo- 
logical defect. Hence, it follows from the above causality argument that during the 
phase transition a network of defects with mean separation £(£) <t will form. 

For models of structure formation and for defects formed in grand unified 
phase transitions we are interested in models with a scale of symmetry breaking 
r\ ~ 10 16 GeV corresponding to a time of formation t c ~ 10 _35 sec. 

The evolution of the field configuration (and thus of the network of defects) 
may be very complicated. However, the same causality argument as used above 
tells us that no correlations on scales > t can be established (see, however, the 
caveats of Refs. 70 and 71). Hence, even at times t >> t c , a network of defects will 
persist with 



For some applications, it is important to have a realistic estimate for the initial 
separation of defects rather than the general causality bound (4.1). In order to 
improve on (4.1) we must assume that ip is in thermal equilibrium. Our goal is to 
estimate £(£) once the defects have stabilized after the phase transition. 



Z(t c ) < t, 



(4.1) 



m < t. 



(4.2) 
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rf)^ll-(^) I • (4-3) 



We follow the methods of Refs. 1,3 and 72. The first step of the analysis is 
to calculate the location <p m (T) of the minimum of the finite temperature effective 
potential Vt (</?)■ From (3.16) it follows that for T <T C 

Tc 

The difference in free energy between (p = (p m (T) and p = is 

AV(T) = V(0) - V T {fp m {T)) = \\^ (l - (|) ' j . (4.4) 

The proto-domain size is given by equating spatial gradient and potential energy 
density 



2 x 2 



giving 

-1/2 

(4.6) 



e(T),A-V V i 1_(0 



The Ginsburg temperature Tq is defined as the temperature below which a thermal 
fluctuation has insufficient energy to take a correlation volume of 'true vacuum' 
V 9 — <Pm{T) into the symmetric configuration ip = 0. Thus, the criterion for Tq is 

e 3 (T G )A\/(T G ) = Tq. (4.7) 

Inserting (4.4) and (4.6) we obtain 

(2\ I/ 2 
(4.8) 

and therefore (from (4.6)) 

i{T G ) = X- 1 ^ 1 ~ A" V 1 . (4.9) 
As a result of the above analysis, it follows that the network of topological defects 
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'freezes out' (i.e. is not continuously destroyed and recreated by thermal fluctua- 
tions) at time t@ corresponding to the temperature Tq. The correlation length at 
the Ginsburg temperature is microscopic, i.e. proportional to i]^ 1 , and thus much 
smaller than the naive causality argument (4.1) would indicate. This fact will be 
important for cosmic string driven baryogenesis. 

Note that the Kibble mechanism was discussed above in the context of a global 
symmetry breaking scenario. As pointed out in Ref. 73, the situation is more 
complicated in local theories in which gauge field can cancel spatial gradients in 
(p in the energy functional, and in which spatial gradients in (p can be gauged 
away. Nevertheless, as demonstrated numerically (in 2 + 1 dimensions) in Ref. 
74 and shown analytically in Ref. 75, the Kibble mechanism also applies to local 
symmetries. 

4.2. Domain Wall and Monopole Problems 

As stated in Table 1, models with domain walls and local monopoles are ruled 
out on cosmological grounds. In both cases, the problem is that the energy density 
in the defects dominates over the energy density in radiation already at the time 
of nucleosynthesis; in other words, the defects would overdose the Universe. The 
reasons why this problem occurs, however, are different for domain walls and local 
monopoles. 

Let us demonstrate explicitly why stable domain walls are a cosmological 
disaster 76 ^ . If domain walls form during a phase transition in the early Universe, it 
follows by causality (see however the caveats of Refs. 70 and 71) that even today 
there will be at least one wall per Hubble volume. Assuming one wall per Hubble 
volume, the energy density pdw °f matter in domain walls is 



PDw(t) ~ V t 



,3 + -l 



(4.10) 



whereas the critical density p c is 



H 



2 



3 




-2 




Pc 



8nG 
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Hence, for r) ~ 10 16 GeV the ratio of (4.10) and (4.11) is 



The above argument depends in an essential way on the dimension of the defect. 
One cosmic string per Hubble volume leads to an energy density p cs in string 

p cs ~ri 2 t~ 2 , (4.13) 

which scales like the background radiation density. Later in this chapter we shall 
see that the scaling (4.13) does indeed hold in the cosmic string model. Hence, 
cosmic strings do not lead to cosmological problems. In the contrary, since for 
GUT models with rj ~ 10 16 GeV 



Pcs - ' V * ~KT 6 , (4.14) 



Pc \m pe 

cosmic strings in these models could provide the seed perturbations responsible for 
structure formation. 

Theories with local monopoles are ruled out on cosmological grounds 77 ) (see 
again the caveats of Refs. 70 and 71) for rather different reasons. Since there are 
no long range forces between local monopoles, their number density in comoving 
coordinates does not decrease. Since their contribution to the energy density scales 
as a~ 3 (t), they will come to dominate the mass of the Universe, provided r\ is 
sufficiently large. 

Theories with global monopoles 78 -' are not ruled out, since there are long range 
forces between monopoles which lead to a "scaling solution" with a fixed number 
of monopoles per Hubble volume. 
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4.3. Cosmic String Evolution 



Applied to cosmic strings, the Kibble mechanism implies that at the time of 
the phase transition, a network of cosmic strings with typical step length will 
form. According to numerical simulations 79 -', about 80% of the initial energy is in 
infinite strings and 20% in closed loops. 

The evolution of the cosmic string network for t > to is complicated . The 
key processes are loop production by intersections of infinite strings (see Fig. 17) 
and loop shrinking by gravitational radiation. These two processes combine to 
create a mechanism by which the infinite string network loses energy (and length 
as measured in comoving coordinates). It will be shown that as a consequence, the 
correlation length of the string network is always proportional to its causality limit 

S(t)~t. (4.15) 

Hence, the energy density Poo(^) i n long strings is a fixed fraction of the background 
energy density p c {t) 

p<x,(f) ~ ^(ty 2 ~ fit- 2 (4.i6) 



or 



(4.17) 

Pc{t) * 1 ; 



Figure 17. Formation of loops by self intersection of infinite strings. According 
to the original cosmic string scenario, loops form with radius R determined by the 
instantaneous correlation length of the infinite string network. 
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We conclude that the cosmic string network approaches a "scaling solution" 80 ) 
in which the statistical properties of the network are time independent if all dis- 
tances are scaled to the horizon distance. 

The origin of the scaling solution for the infinite string network can be under- 
stood heuristically as follows. If the curvature radius £ (t) of this network is much 
larger than the Hubble radius t, the network will be frozen in comoving coordinates 
(since the Hubble damping term dominates in the equations of motion). Hence in 
the radiation dominated FRW epoch 

f(t) ~ a(t) ~ t 1 / 2 (4.18) 

and the Hubble radius will catch up to £(£). Conversely, if <C t then the tension 
term in the equations of motion for the string will dominate, the strings will oscillate 
relativistically and there will be frequent self intersections of the strings, leading to 
rapid loop production and to increasing £(t)/t. Combining these two arguments, 
we conclude that there must be a 'dynamical fixed point' with £(t) ~ t. 

A first step in a more rigorous analysis of cosmic string evolution is the deriva- 
tion of the effective equation of motion for the strings, Note that this equation must 
follow from the field equations since the string is merely a particular topologically 
stable field configuration. 

The equations of motion of a string can be derived from the Nambu action 

/l /2 
dadr (- det g^j a, b = 0,1, (4.19) 

(2) 

where g ab is the world sheet metric and a and r are the world sheet coordinates. In 
flat space-time, r can be taken to be coordinate time, and a is an affine parameter 
along the string. In terms of the string coordinates X^(a,r) and the metric gffi 
of the background space-time, 

9$ = X^JtS • (4.20) 
From general symmetry considerations, it is possible to argue that the Nambu 



45 



action is the correct action. However, I shall follow Foerster 81 ) and Turok 82 ) and 
give a direct heuristic derivation. We start from a general quantum field theory 
Lagrangian Cqft- The action is 




(4.21) 



We assume the existence of a linear topological defect at X IJ, (a, r). The idea 
is to change variables so that o and r are two of the new coordinates, and to 
expand S to lowest order in w/R, where w is the width of the string and R its 
curvature radius. As the other new coordinates we take coordinates p 2 and p 3 
in the normal plane to X fi (a,r). Thus the coordinate transformation takes the 
coordinates y M (/i = 0, 1, 2, 3) to new ones a v = (r, a, p 2 , p a ): 



where i — 2, 3 and are the basis vectors in the normal plane to the string world 
sheet. The measure transforms as 



y 



^a a )=X^a,T)+p*nf(a,T) 



(4.22) 




(4.23) 



with 




(4.24) 



The determinant can easily be evaluated using the following trick 



detM£ = (-det^M^M^) 172 = ^-det D, 



ab 



(4.25) 




(4.26) 



4G 



Hence 




(4.27) 



Here, —fj, is the integral of C in the normal plane of X. To first order in w/R, it 
equals the integral of —7i; hence it is the mass per unit length. 

This derivation of the Nambu action is instructive as it indicates a method for 
calculating corrections to the equations of motion of the string when extra fields are 
present, e.g. for superconducting cosmic strings. It also gives a way of calculating 
the finite thickness corrections to the equations of motion which will be important 
at cusps (see below). 

In flat space-time we can consistently choose r = t, x ■ x' — and x 2 + x' = 0. 
The equations of motion derived from the Nambu action then become 



where / indicates the derivative with respect to a. As expected, we obtain the 
relativistic wave equation. The general solution can be decomposed into a left 
moving and a right moving mode 




(4.28) 




(4.29) 



The gauge conditions imply 




(4.30) 



For a loop, x(<7, t) is periodic and hence the time average of a and b vanish. Thus, 
a and b are closed curves on the unit sphere with vanishing average. Two such 
curves generically intersect if they are continuous 83 -* . An intersection corresponds 
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to a point with x' = and x = 1. Such a point moving at the speed of light is 
called a cusp. Note that x(a, t) need not be continuous. Points of discontinuity 
are called kinks. Both cusps and kinks will be smoothed out by finite thickness 
effects 84 ). 

The Nambu action does not describe what happens when two strings hit. This 
process has been studied numerically for both global 85 ) and local 86 ) strings. The au- 
thors of those papers set up scalar field configurations corresponding to two strings 
approaching one another and evolve the complete classical scalar field equations. 
The result of the analysis is that strings do not cross but exchange ends, provided 
the relative velocity is smaller than 0.9. Thus, by self intersecting, an infinite string 
will split off a loop (Figure 17). An important open problem is to understand this 
process analytically. For a special value of the coupling constant Ruback has given 
a mathematical explanation 87 ) (see also Shellard and Ruback in Ref. 86). 

There are two parts to the nontrivial evolution of the cosmic string network. 
Firstly, loops are produced by self intersections of infinite strings. Loops oscillate 
due to the tension and slowly decay by emitting gravitational radiation 88 ). Com- 
bining the two steps we have a process by which energy is transferred from the 
cosmic string network to radiation. 

There are analytical indications that a stable "scaling solution" for the cosmic 
string network exists 4 ) . It is given by on the order 1 infinite string segment crossing 
every Hubble volume. The correlation length £(£) of an infinite string is of the order 
t. Hence, at time t loops of radius R ~ t are produced, of the order 1 loop per 
Hubble volume per expansion time. A heuristic argument for the scaling solution 
is due to Vilenkin. Take v(t) to be the mean number of infinite string segments 
per Hubble volume. Then the energy density in infinite strings is 

Poo (t) =fii)(t)t- 2 (4.31) 

The number of loops n(t) produced per unit volume is given by 

= CZ>V4 (432) 
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where c is a constant of the order 1. Conservation of energy in strings gives 

^PoqQ) 3 , dn , _ 2 3 

~dT + 2t = ~^t^ = -c^t (4.33) 

or, written as an equation for v{t) 

I - ^ = -cc'vh- 1 (4.34) 

Thus if v ^> 1 then z> < while if v <C 1 then z> > 0. Hence there will be a stable 
solution with i> ~ 1. 

The precise value of i> must be determined in numerical simulations. These 
simulations are rather difficult because of the large dynamic range required and due 
to singularities which arise in the evolution equations near cusps. In the radiation 
dominated epoch, v is still uncertain by a factor of about 10. The first results 
were reported in Ref. 89. More recent results are due three groups. Bennett and 
Bouchet 90 -* and Allen and Shellard 91 ) are converging on a value 10 < v < 20, 
whereas Albrecht and Turok 92 ) obtain a value which is about 100. 

4.4. Scaling Solution for Strings 

The scaling solution for the infinite strings implies that the network of strings 
looks the same at all times when scaled to the Hubble radius. This should also im- 
ply that the distribution of cosmic string loops is scale invariant in the same sense. 
At present, however, there is no convincing evidence from numerical simulations 
that this is really the case. 

A scaling solution for loops implies that the distribution of Ri(t), the radius of 
loops at the time of formation, is time independent of time after dividing by t. To 
simplify the discussion, I shall assume that the distribution in monochromatic, i. e. 

Ri(t)/t = a. (4.35) 

From Figure 17, we expect a ~ 1. The numerical simulations, however, now give 
a < 10- 2 90 ' 91 ). 
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From the scaling solution (4.15) for the infinite strings we can derive the scaling 
solution for loops. We assume that the energy density in long strings - inasmuch 
as it is not redshifted - must go into loops. (5 shall be a measure for the mean 
length I in a loop of "radius" R 

e = (3R. (4.36) 

Since per expansion time and Hubble volume about 1 loop of radius Ri(t) is pro- 
duced, we know that the number density in physical coordinates of loops of radius 
Ri(t) is 

n(Ri(t),t) = ct- 4 (4.37) 

with a constant c which can be calculated from (4.31), (4.35) and (4.36). Neglecting 
gravitational radiation, this number density simply redshifts 

n(R,t)= (^^y) n(R,t f (R)), (4.38) 

where t f(R) is the time when loops of radius R are formed. Isolating the R 
dependence, we obtain 

n(R,t) ~ R~ 4 z(R)~ 3 (4.39) 

where z(R) is the redshift at time t = R. We have the following special cases: 

n(R, t) ~ iT 5/ V 3/2 t<t eq 
n(R, t) ~ R-V 2 t l J q 2 t- 2 t>t eq ,t f (R)< t eq (4.40) 
n(R, t) ~ iT 2 t~ 2 t>t eq ,tf(R)> t eq . 

The proportionality constant c is 

c=^p- 1 a- 2 i> (4.41) 

(see e.g. Ref. 93). In deriving (4.41) it is important to note that n(Ri(t),t)dRi 
is the number density of loops in the radius interval [Ri, Ri + dRi}. Hence, in the 
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radiation dominated epoch 

n(R, t) = vR-Vh-V 2 



(4.42) 



with 

u = ]-(3' l a l/2 u . (4.43) 

From (4.43) we can read off the uncertainties in v based on the uncertainties 
in the numerical results. Both a 1 / 2 and v are determined only up to one order of 
magnitude. Hence, any quantitative results which depend on the exact value of v 
must be treated with a grain of salt. 

Gravitational radiation leads to a lower cutoff in n(R,t). Loops with radius 
smaller than this cutoff were all formed at essentially the same time and hence have 
the same number density. Thus, n(R) becomes flat. The power in gravitational 
radiation Pq can be estimated using the quadrupole formula 94 ). For a loop of 
radius R and mass M 

P G = \g<QQ>, (4.44) 

where Q is the quadrupole moment, Q ~ MR 2 , and since the frequency of oscilla- 
tion is to = R^ 1 

P G ~ G(MR 2 ) 2 uj 6 ~ (G/i)/i . (4.45) 

Even though the quadrupole approximation breaks down since the loops move 
relativistically, (4.45) gives a good order of magnitude of the power of gravitational 
radiation. Improved calculations 88 -* give 

P G = 7 (G/i)/i (4.46) 

with 7 ~ 50. (4.36) and (4.46) imply that 

R = (4.47) 

with 7 = 7//? ~ 5 (using /3 ~ 10 89 )). Note that the rate of decrease is constant. 
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Hence, 



R(t) =Ri-(t- ttfGn 



(4.48) 



and the cutoff loop radius is 



R c ~ jGfiU . 



(4.49) 



Let us brieffy summarize the scaling solution: 

f ) At all times the network of infinite strings looks the same when scaled by the 
Hubble radius. A small number of infinite string segments cross each Hubble 
volume and poo(t) is given by (4.31). 

2) There is a distribution of loops of all sizes < R < t. Assuming scaling for 
loops, then 



Although the qualitative characteristics of the cosmic string scaling solution 
are well established, the quantitative details are not. The main reason for this 
is the fact that the Nambu action breaks down at kinks and cusps. However, 
kinks and cusps inevitably form and are responsible for the small scale structure 
on strings. In fact, coarse graining by integrating out the small scale structure 
may give an equation of state for strings which deviates from that of a Nambu 
string 95 ^. Attempts at understanding the small scale structure on strings are at 
present under way 96 ). 




(4.50) 



where a 



R is the time of formation of a loop of radius R. Also 



n(R, t) = ntfGfrt, t) , R < jGfxt . 



(4.51) 
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4.5. Scaling Solution for Textures 

Applied to textures 69 ), the Kibble mechanism implies that on all scales r > 
h e ld configurations with nonvanishing winding number frozen in after 

the phase transition at time tQ. 

The general causality argument implies that at all times t > to, on scales 
greater than t the field configuration will be random, and hence configurations 
with winding number n w ^ will persist. 

The texture dynamics is determined by the field equation 

<p + 3Hip-a- 2 V 2 (p = -V'(if). (4.52) 

The second term on the left hand side is proportional to H 2 tp, the third term 
scales as k 2 a~ 2 p, where k' 1 is the wavelength of the inhomogeneity of the field 
configuration. Comparing these two terms, we see that for 

ha" 1 < H, (4.53) 

i.e. on scale larger than the Hubble radius, a texture configuration will be frozen 
in, whereas for 

ka' 1 > H, (4.54) 

on scales smaller that the Hubble radius, the motion of ip will be relativistic, and 
the texture configuration will become homogeneous on a time scale of t. 

Hence, in the texture scenario, at time t ^> to, the field configuration tp(x) 
is predicted to be homogeneous on scales smaller than t, but inhomogeneous on 
larger scales. There is a finite (and time independent) probability p(n c ) that per 
horizon volume there is a field configuration with winding n w greater than n c . The 
probability p(n c ) can be determined using combinatorial arguments 97 \ 
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The evolution of texture configurations has been studied numerically in Refs. 
98-100 and analytically in Ref. 101 (see also Ref. 102). The qualitative features of 
the evolution depend on wether n w is greater or smaller than the critical winding 
n c . If n w > n c , then the texture will contract, the winding will increase, and 
eventually the texture will unwind by (p jumping over the potential barrier at the 
texture center. On the other hand, if n w < n c , the field configuration will dissipate. 
It will expand with n w decreasing to zero continuously. 

The physics of texture evolution can be most readily analyzed by studying the 
spherically symmetric configuration (3.55) with x a function of r and t. For such 
a configuration, the action in an expanding FRW background metric is 

S = J d 4 x r 2 a 3 [x 2 - a- 2 X ' 2 - 2a~ 2 V~ 2 sin 2 x ]- (4.55) 

The second term corresponds to the radial gradient energy, the third to angular 
gradient energy. 

Consider a field configuration as sketched in Fig. 18 a with Xmax > n c > ir/2. 
In this case, the forces F\ and F<i act as indicated in the figure. In order to 
minimize the angular gradient energy, Xmax wants to increase (see Fl). However, 
for r < x there is an additional (and in part opposing) force F% in order to reduce 
the angular gradient energy, r(x) will tend to decrease, i.e. the field configuration 
will contract and its total winding will increase. 

However, for Xmax < 7r/2, the angular gradient force tends to reduce Xmax, 
and to reduce the radial gradients, x( r ) will tend to decrease even for x < Xmax- 
Hence, the field configuration will dissipate (see Figure 18 b). 

If n w is only slightly larger than 7r/2, the force Fl of Fig. 15 a is too weak to 
offset the force F2 described above. The critical winding n c is hence larger than 
0.5. The precise value has been determined in Refs. 99 and 101. It depends on the 
length of the 'plateau region' of x( r ) (which is the inter-texture separation) and 
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on the expansion rate of the Universe. For realistic parameters 101 ) 

0.65 <n c < 0.75 . (4.56) 

As will be shown in a later Chapter, only textures with n w > n c generate localized 
density perturbations which can act as seeds for cosmic structure formation. 



Figure 18. A sketch of the forces acting on a radially symmetric texture configura- 
tion and which cause unwinding in Fig. 18a if n w > n c , and dissipation if n w < n c 
(Fig. 18b). 



5. Introduction to Structure Formation 

5.1. Power Spectrum 

In Chapter 2, the cosmology of a homogeneous and isotropic Universe was 
reviewed. In order to understand structure formation, it is essential to study 
the evolution of inhomogeneities at a linearized level. This will be adequate to 
understand the early evolution of density perturbations in the Universe. 

Starting point of the relativistic theory of cosmological perturbations 103 ) are 
the linearized Einstein equations. If we take the general Einstein equations 

GV = SttG 7>, (5.1) 

where G^ v (g a p) is the Einstein tensor and T^ v is the energy-momentum tensor of 
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matter, and expand about a cosmological background solution 

g^ = diag (l, -a 2 (t), -a 2 (t), -a 2 (t)) (5.2) 

and 

TP® = diag(p,- P ,- P ,-p) (5.3) 
(see Chapter 2), then we obtain the linearized equations 

SG^ (<? Q /°) + h a(3 ) = 8nG ST^. (5.4) 

Equation (5.4) relates the perturbation h^ v of the metric, i.e. 

to the matter perturbations. 

To gain a heuristic understanding of how the perturbations evolve, recall that 
gravity is a purely attractive force. Given an initial mass perturbation 5m, the 
force on surrounding particles will be 

F ~ 5m (5.6). 

Since (neglecting for a moment the expansion of the Universe) 

5m ~ F, (5.7) 

we see that in a nonexpanding background the growth of perturbations is expo- 
nential. In an expanding background, there will be a damping term depending on 
H . Hence, perturbations will increase only as a power of time. 
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The details of the analysis are rather complicated (see Refs. 30 and 104 for 
recent reviews). The result is that the density contrast 5p grows as follows 

f t 2/3 t > t eq 
5p(t) ~ I t t<t eq , \ >t (5.8) 
I, const t < t eq , A < t 

Note that on length scales A greater than the Hubble radius t, the quantity bp is 
not gauge invariant, i.e. it depends on the slicing of space-time 104 ). The quantity 
which is gauge invariant is the relativistic potential $, which is time independent 
if the equation of state of the background cosmology is constant 104 ). In a gauge 
in which g^y is diagonal and for models of matter in which 5T{j is diagonal at 
linearized level (a condition satisfied by most interesting models of matter), <3> can 
be identified as follows: 

g^ = (1 + 2$)dt 2 - a 2 (t)(l - 2<5>)dx 2 . (5.9) 

We will use the results of (5.8) when describing the evolution of the power spectrum. 

The main quantity of interest is the r.m.s. mass excess (SM/M) (k,t) at time 
t in a sphere of radius fc -1 . Given a smooth density distribution 

p(x,t)=po(t) + 6p(x,t), (5.10) 

the r.m.s. mass excess can be related to the Fourier mode Sp(k) in a straightforward 
manner 18 -*. The result is 



SM\ n _ ^ _ 7 3 1 5p , 2 



2 

(M)^ 3 Ff(M). (5.H) 
M J po 

The adopted convention for Fourier transformation is 

5p(x) = (2tt)- 3 /V- 1 / 2 j d 3 ke^5p(k) . (5.12) 

The result (5.11) holds provided \5p(k)\ 2 is proportional to k n with n > —3. An 
intuitive way to understand the result is as follows: perturbations with wave num- 
ber larger than k average to zero in a volume A; -3 , perturbations with wave number 
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smaller then k are phase space suppressed such that (5M/M) (k) receives its major 
contribution from Fourier modes of wave number k. Their phase space volume is 
k 3 . 

The most commonly used function describing the ensemble of perturbations is 
the spectrum P(k). By definition, the power spectrum is the square of the modulus 
of the Fourier space density contrast 

P{k) = \ 6 Ak)\\ (5.13) 
Po 

Hence (from (5.11)), P(k) is related to the r.m.s. mass fluctuations (SM/M) (k,t) 
on physical length scale 

\k = a(t) T (5.14) 

at time t via 

(^fj (k,t)^k*P(k). (5.15) 

The scaling solution for topological defect models implies that when measured 
at the time tu(k) when the wavelength equals the Hubble radius, the r.m.s. 
mass perturbation (SM/M) (k) is independent of k, i.e., 

^L(k,t H (k)) = const. (5.16) 

This is because at any time t, a constant fraction of the mass M inside the Hubble 
radius is contained in the topological defects. For example, one cosmic string of 
length t contains mass SM = fit compared to the total mass M ~ t 3 p(t) ~ t inside 
the Hubble radius, thus leaving the ratio SM/M time independent. 

Eq. (5.16) is the same result as is obtained for inflationary Universe models. 
Hence, we conclude that all three main models of structure formation: adiabatic 
random phase perturbations from inflation, cosmic strings, and global textures, to 
a first approximation produce a scale-invariant spectrum. 
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To convert (5.16) into an expression for the power spectrum P(k), we use the 
fact that 5M/M grows as the scale factor a(t) during the matter dominated epoch 
on scales smaller than the Hubble radius (see (5.8)) 



5M . . ft \ 2/3 SM 



m (m > = UmJ ■*<*•*<*»■ (5 - 17) 

On scales larger than the Hubble radius at t eq 

t H {k) = 2nk- 1 a(t H (k)) ~ t 2 J?(k) k~ x , (5.18) 

and hence 

t H (k)~k~ 3 . (5.19) 
Therefore, combining (5.16), (5.17) and (5.19) 

8M 

From (5.15) it follows that 



(k,t)~k 2 . (5.20) 



P(k) ~ fc n (5.21) 



with n — 1. 



Recently, there has been some interest in deviations from scale invariance. In 
models of inflation, a deviation comes about 105 ) because H decreases slowly during 
inflation. In topological defect models, numerical 106 ' 107 -* and semi-analytical 108 -* 
studies have also shown small deviations from scale invariance. These deviations, 
however, all small and quite model dependent. 
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5.2. CMB Anisotropies 



Density perturbations give rise to anisotropies in the temperature of the CMB. 
There are three main contributions (see Fig. 19): 

i) Density fluctuations at the time of last scattering t rec imply that the sur- 
face of constant temperature T rec was not flat. Hence, photons arriving at 
the observer O from different directions (separated by angle 9) have trav- 
elled a different amount of time. Since they were emitted with the same 
frequency, they arrive with differing frequencies, thus leading to temperature 
anisotropies. This is the Sachs- Wolfe effect 109 ). 

ii) A localized clump of energy between last scattering and the present time 
will distort the geodesies passing through it, and hence lead to temperature 
differences (the Rees-Sciama effect 110 )). 

Hi) If either the last scattering surface or the observer are moving with respect 
to the frame given by the background cosmology, Doppler anisotropies will 
result. 



Figure 19. Space-time plot sketching the origin of CMB temperature anisotropies. 
The surface labelled T rec is the last scattering surface. O is the observer at the 
present time measuring photons 7 impinging from directions in the sky separated 
by angle The shaded area labelled C stands for a local overdensity, leading to 
distortions of geodesies. Possible velocities of observer and emitter are indicated as 
v and v e . 
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The observed m ) dipole anisotropy of the CMB is used to determine the 
Earth's velocity relative to the rest frame of the CMB. After subtracting the dipole, 
the remaining anisotropy on large angular scales is dominated by the Sachs- Wolfe 
effect (see e.g. Ref. 30). The specific relationship between 5T/T and density 
perturbations in the case of adiabatic fluctuations is given by 30 > 104 ) 

^(f o ,0) = £*(f rec ,A(0)). (5.22) 

According to relativistic perturbation theory, the relativistic potential $ of (5.9) 
is independent of time in the absence of changes in the equation of state of the 
background 104 -'. Hence 



— (i ,tf) = 3*(te(A),A(0)) ~ -^(t H (X),X) (5.23) 

where A(i?) is the co moving length scale corresponding to angular separation 6 at 
t re c, and iff (A) is the time when this scale equals the Hubble radius i. We have 
used the fact that at Hubble radius crossing <3> and Sp/p coincide up to a numerical 
factor of the order 1 . 104 ) 

Recently, both the COBE satellite n ) and MIT balloon 112 ) experiment have 
detected CMB anisotropies on angular scales larger than 7°. The quadrupole 
anisotropy reported is 

^(90°) ~ 5 • 10" 6 , (5.24) 

and the low / harmonies of 5T/T are consistent with a power law spectrum of 
density perturbations with 

n= 1.1 ±0.5. (5.25) 

It is reassuring that the observed value of n is so close to the value n = 1 favored 
in both inflationary and topological defect models of structure formation. The 
amplitude of 5T/T is also of the right order of magnitude to agree with theoretical 
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predictions based on structure formation arguments. Since by (5.8) density per- 
turbations grow as a(t) and since today Sp/p on scales of clusters is about 1, we 
predict that on these scales 

^W~^(W _1 ~3.10- 5 . (5.26) 

This result agrees with the extrapolation of the COBE results to cluster scales n ) 
to better than a factor of 2. 

The results (5.24) and (5.25) cannot be used to differentiate between inflation- 
ary Universe and topological defect theories, since the models all predict a similar 
slope of P(k). However, the amplitude of 5T/T can be used to normalize the 
power spectrum in any given model. Large scale structure observations provide an 
independent normalization. These two normalizations must be consistent in order 
for a theoretical model to work. For cosmic strings, the normalizations of P(k) 
agree 106 ' 108 -* well, for an unbiased texture model, the normalization factors differ 
by about 3 107 ). 

In the near future, maps of CMB anisotropies will be obtained 113 ) which are 
signal dominated in every pixel (the COBE maps are dominated by noise). At that 
point, statistical 114 ' 115 ^ measures of CMB maps can be evaluated which can pick 
out the non-Gaussian signatures of topological defect models. Indeed, non-random 
phases of the Fourier modes of 5p are the key feature of such models (see Fig. 20). 

Note that CMB anisotropies form a more severe testing ground for cosmological 
models than large-scale structure data obtained from optical telescopes because 
they are (at least an angular scales larger than the horizon at last scattering) not 
contaminated by hydrodynamical and nonlinear effects. 
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Figure 20. A comparison of density perturbations along a line in space between 
(a) random phase models and (b) theories based on nonadiabatic seeds. 

5.3. Large- Scale Structure Data 

Here, I shall briefly mention some of the data obtained from optical telescopes 
which at the present time seems most relevant to cosmology. 

First, however, it is important to note a discrepancy between theory and obser- 
vations. Observations become increasingly uncertain on increasing length scales, 
whereas theoretical predictions decrease in accuracy as the length scale decreases, 
the reason being that on large scales linear theory is applicable and gravity is the 
only force which needs to be considered, whereas on smaller scales nonlinear and 
messy nongravitational effects become important. However, this state of affairs 
also implies that progress either on the theoretical or observational front will yield 
double benefit. For example, new data on larger scales not only increases the 
amount of data available, but also leads into a regime with smaller theoretical un- 
certainties. Since we know that observers will be providing lots of new data in the 
next few years, we can be sure that cosmology will remain an extremely exciting 
field. 

There is a lot of new data on the large-scale structure of the Universe. Red- 
shift surveys have provided three dimensional maps of the distribution of galax- 
ies. An example is the recent Center for Astrophysics (CFA) survey 9 ) of slices 
on the northern celestial sphere which show overdense sheets of galaxies with di- 
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mensions (50 x 50 x 5)Mpc (with h — 1) separated by large voids. There is 
evidence for superclusters 116 -* , filaments longer than 50/i _1 Mpc 117 ), voids of di- 
ameter 60/i _1 Mpc 118 ^ and overdense sheets of galaxies 119 -* in regions of the sky 
different from that covered in the CFA survey. 

A second window on the large-scale structure of the Universe comes from 
measuring the peculiar velocities of galaxies averaged over large regions 120 -* from 
which one can infer information about the magnitude of density perturbations on 
large scales 121 -*. Preliminary observations indicate surprisingly large velocities on 
scales of 60/i _1 Mpc, although consensus on the interpretation of the data is still 
lacking 120 - 122 ). 

On the scale of clusters the data is still rather uncertain. The quantities one 
would like to focus on is the overall mass scale (the mean mass of an object which 
satisfies a fixed operational definition of a cluster), the multiplicity function and 
the two point (and in the future higher point) correlation functions. 

A cluster of galaxies (Abell 123 ) cluster) is defined as a region in the sky with 
more than 50 bright galaxies in a sphere of radius 1.5/i _1 Mpc. (To be compared 
to the mean separation of bright galaxies 124 ) which is about 5/i _1 Mpc). The 
mean separation of clusters is about 50/i _1 Mpc, and their masses are of the order 
1O 14 M0 125 ). 

The multiplicity function n(M) of clusters gives the number density of clusters 
of mass M per unit physical volume. n(M)dM is the number density of objects in 
the mass interval [M, M + dM]. 

The two point correlation function £(r) of clusters measures the nonrandomness 
of the distribution of clusters and is defined by 

^) = <IiM^£>, (5.27) 
no 

where n(r) is the number density at a fixed distance r from a given object and no 
is the average density, and pointed brackets indicate averaging over the objects. 
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The initial observational results 10 ) (Figure 24) gave 

, N-1.8 

i{r) ~ ( (5.28) 

with correlation length tq ~ 25/i _1 Mpc. Recently 126 ) , some criticism of these 
results has been raised. However, it is unlikely that the entire effect is fictitious. 



Figure 21. The observed cluster correlation function plotted against relative sep- 
aration (from Ref. 10). The solid circles are the data points. Shown is the fit to a 
r~ 2 slope. 

For galaxies, similar quantities can be measured and compared with observa- 
tions. Large spiral galaxies have masses M ~ 1O U M0. Their mean separation 
is about 5/i _1 Mpc. The mass function 127 ) (Figure 22) gives the distribution of 
masses, and the two point correlation function describes the nonrandomness of the 
distribution. The two point correlation function of galaxies has the same power 
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Figure 22. The mass function of galaxies (determined from the luminosity func- 
tion assuming constant mass to light ratio) from Refs. 127 ((a) from Bahcall, (b) 
from Binggeli). 

law form as (5.28), with correlation length tq ~ 8Mpc. For galaxies, one can in 
addition probe the interior structure by measuring the velocity rotation curve v(r) 
which is related to the density profile (baryonic plus dark matter) by 

v(r) ~ r 2 p(r) . (5.29) 

Observations 128 ^ indicate that far beyond the disk radius, v(r) ~ const (Figure 23) 
which indicates the existence of dark matter with a density profile 

p(r)~r~ 2 . (5.30) 

We can also measure the angular momentum of galaxies. Typical numbers for large 
spirals are in the range 129 ) 10 73 — 10 75 cm?g s _1 . 
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Figure 23. A typical velocity rotation curve (for NGC488 HI data, taken from 
Ref. 128). The radius is in kpc, the velocity in kms~ l . 

6. Cosmic Strings and Structure Formation 

Starting point of the structure formation scenario in the cosmic string theory is 
the scaling solution for the cosmic string network, according to which at all times 
t (in particular at t eq , the time when perturbations can start to grow) there will 
be a few long strings crossing each Hubble volume, plus a distribution of loops of 
radius R < t (see Fig. 24). 



Figure 24. Sketch of the scaling solution for the cosmic string network. The box 
corresponds to one Hubble volume at arbitrary time t. 
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The cosmic string model admits three mechanisms for structure formation: 
loops, filaments and wakes. Cosmic string loops have the same time averaged field 
as a point source with mass 130 ) 

M{R)=(3Rn, (6.1) 

R being the loop radius and (3 ~ 2n. Hence, loops will be seeds for spherical 
accretion of dust and radiation. 80 ) 

For loops with R <t eq , growth of perturbations in a model dominated by cold 
dark matter starts at t eq . Hence, the mass at the present time will be 

M(R,t ) = z(t eq )PRfi. (6.2) 

In the original cosmic string model 80,131 ) it was assumed that loops dominate 
over wakes. In this case, the theory could be normalized (i.e., /i could be deter- 
mined) by demanding that loops with the mean separation of clusters d c \ (from 
the discussion in Section 4.4 it follows that the loop radius R(d c i) is determined 
by the mean separation) accrete the correct mass, i.e., that 

M(R(d d ),t ) = 10 u M Q . (6.3) 

This condition yields 131 ) 

fi ~ 10 36 GeV 2 (6.4) 

Thus, if cosmic strings are to be relevant for structure formation, they must arise 
due to symmetry breaking at energy scale r\ ~ 10 16 GeV. This scale happens to 
be the scale of unification of weak, strong and electromagnetic interactions. It is 
tantalizing to speculate that cosmology is telling us that there indeed was new 
physics at the GUT scale. 
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Figure 25. Sketch of the mechanism by which a long straight cosmic string moving 
with velocity v in transverse direction through a plasma induces a velocity perturba- 
tions Av towards the wake. Shown on the left is the deficit angle, in the center is a 
sketch of the string moving in the plasma, and on the right is the sketch of how the 
plasma moves in the frame in which the string is at rest. 

The second mechanism involves long strings moving with relativistic speed in 
their normal plane which give rise to velocity perturbations in their wake. The 
mechanism is illustrated in Fig. 25: space normal to the string is a cone with 
deficit angle 132 ) 

a = 8nG/i. (6.5) 

If the string is moving with normal velocity v through a bath of dark matter, a 
velocity perturbation 

5v = 47rG/it>7 (6.6) 

[with 7 = (1 — t> 2 )~ 1//2 ] towards the plane behind the string results 133 ) . At times 
after t eq , this induces planar overdensities, the most prominent (i.e. thickest at the 
present time) and numerous of which were created at t eq , the time of equal mat- 
ter and radiation 134 ' 135 ^. The corresponding planar dimensions are (in comoving 
coordinates) 

t eq z(t eq ) x t eq z(t eq )v ~ (40 x 40*;) Mpc 2 . (6.7) 

An intuitive understanding of the origin of the above distinguished scale can 
be obtained as follows. Viewed from a distance, the density perturbation grows as 
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in a linear theory i.e. 

^( t) = (i) 2/ %) (6.8) 
P \UJ P 

for a perturbation set up at a time ti > t eq . Since the initial 5p/ p(ti) is independent 
of ti, the largest density contrast comes from the earliest ti, namely U = t eq . 
Fluctuations created at U < t eq are erased by the large thermal velocities. Thus, 
strings at t eq create the most prominent wakes. They are also the must numerous, 
since the comoving separation of strings decreases as ti decreases. 

A more rigorous way to obtain the above result is to consider the evolution of 
the velocity perturbations induced by a wake in the Zel'drovich approximation 30 ). 
The height of a dark matter particle above the wake can be written as 

h(q,t) = a(t)[q-i/>(q,t)] , (6.9) 

where q is the initial comoving distance, and ip(q,t) is the comoving displacement 
caused by the presence of the wake. The thickness of the wake at time t is deter- 
mined by the value of q for which 

h(q,t) = 0. (6.10) 

Obviously, the value of q for which h — increases as the time when the pertur- 
bation ip begins to grow is moved back in time. Hence, the earliest wakes will be 
the thickest. 

The details of the calculation depend on whether the dark matter is hot or 
cold (see Chapter 8). For hot dark matter, the large thermal velocities delay the 
beginning of the growth of ip on small scales. A detailed analysis 136 > 137 ) shows 
that for hot dark matter no perturbations become nonlinear unless 

Gp>5-10~ 7 (6.11) 

In this case, the value of the redshift z(q) at which h = is maximal for the value 
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of q (the thickness) given by 



d ~ Gfiv 7(u) z(t eq ) 2 t eq ~ 4w Mpc (6.12) 

for wakes created at t{ = t eq . Note that the scales of cosmic string wakes (see (6.7) 

and (6.12)) compare favorably with the measures of the observed sheets of galaxies 
9)_ 

Wakes arise if there is little small scale structure on the string. In this case, 
the string tension equals the mass density, the string moves at relativistic speeds, 
and there is no local gravitational attraction towards the string. 

In contrast, if there is small scale structure on strings, then 95 ) the string tension 
T is smaller than the mass per unit length /x and the metric of a string in z direction 
becomes 

ds 2 = (1 + h 00 )(dt 2 - dz 2 - dr 2 - (1 - 8Gfi)r 2 d0 2 ) (6.13) 

with 

h Q0 = 4G(n-T)]n-, (6.14) 

tq being the string width. Since hoo does not vanish, there is a gravitational 
force towards the string which gives rise to cylindrical accretion, thus producing 
filaments. 

As is evident from the last term in the metric (6.13), space perpendicular to 
the string remains conical, with deficit angle given by (6.5). However, since the 
string is no longer relativistic, the transverse velocities v of the string network are 
expected to be smaller, and hence the induced wakes will be shorter and thinner. 

Which of the mechanisms - filaments or wakes - dominates is determined by 
the competition between the velocity induced by hoo and the velocity perturbation 
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of the wake. The total velocity is 



138,139) 



U = 



2nG(fi - T) 



(6.15) 



the first term giving filaments, the second producing wakes. Hence, for small v the 
former will dominate, for large v the latter. 

By the same argument as for wakes, the most numerous and prominent fila- 
ments will have the distinguished scale 



Starting point of the texture scenario of structure formation is the scaling 
solution for textures: at any time t, there is a fixed probability p that a texture 
configuration is entering the Hubble radius and starting to collapse. 

In the texture model it is the contraction of the field configuration which leads 
to density perturbations. At the time when the texture enters the horizon, an 
isocurvature perturbation is established: the energy density in the scalar field is 
compensated by a deficit in radiation. However, the contraction of the scalar field 
configuration leads to a clumping of gradient and kinetic energy at the center of 
the texture 140 ) (Fig. 26). This, in turn, provides the seed perturbations which 
cause dark matter and radiation to collapse in a spherical manner 140-143 ). 

The texture model has a significant advantage over the cosmic string theory: 
it is much more amenable to analytical and numerical calculations. The reason is 
simple: whereas the evolution of the network of cosmic stings is very complicated, 
and the detailed statistics of the scaling solution are not yet known with good 
accuracy, the evolution of textures is essentially trivial: they collapse in less than 



teqZ{t eq ) X df X df 



(6.16) 



where df can be calculated using the Zel'dovich approximation. 



7. Textures and Structure Formation 
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Figure 26: A sketch of the density perturbation produced by a collapsing texture. 
The left graph shows the time evolution of the field x{ r ) as a function of radius r 
and time (see (3.55)). The contraction of x{ r ) leads to a spatial gradient energy 
perturbation at the center of the texture, as illustrated on the right. The energy is 
denoted by p. Solid lines denote the initial time, dashed lines are at time t + At, 
and dotted lines correspond to time t + 2At, where At is a fraction of the Hubble 
expansion time (the typical time scale for the dynamics). 

an Hubble expansion time and set up well defined cosmological perturbations whose 
subsequent evolution can be analyzed without great problems. 

In the scaling solution, the comoving number density of textures n obeys 8,69 ) . 
where c ~ 0.04 (see also Refs. 97 and 100), and r is conformal time determined by 



dt 



a(t)dr . 



(7.2) 



A collapsing texture with center at x 



induces a velocity perturbation' 



8,140) 



Av = 



F(x, t)ex 



(7.3) 



where x is the unit vector in direction x 



£ = 8vr 2 Gr/ 2 



(7.4) 
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is the measure of the strength of the velocity kick (recall that i] is the symmetry 
breaking scale, the single free parameter in the texture model) and, 

F(x,t)^{ {1 ~^ X<T (7.5) 

I X > T 

gives the profile of the velocity field. In the above x is a comoving coordinate. 
In agreement with the causality conditions, F vanishes on scales larger than the 
horizon (\x\ = x > r). 

The density fluctuation 5(x, r, Tj) at time r induced by a texture collapsing at 
x = at time Ti is determined by the initial velocity perturbation. The result is 

8,141) 

S(x : T,n) = 2F ^i)en a{m)[52{r] . )5l{ri) _ ^(thMt,)] (7.6) 

where Si(r]) and 62 (r/) are the growing and decaying mode solutions of the linear 
perturbation equations, and 



1 _ f8nGp(t eq )\ 1/2 



* 



(7.7) 



In the above, the scale factor has been normalized such that 

a(t eq ) = 1 . (7.8) 

In this case the growing mode is 

5Mv)) = 1 + y ■ (7-9) 

Based on (7.1) and (7.6), many interesting observables can be computed rather 
easily. As an example, in Fig. 27 the resulting power spectrum of density pertur- 
bations is shown and compared to what is obtained in the "standard CDM model" , 
a theory based on quantum fluctuations from inflation in a Universe dominated by 
cold dark matter. 
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Figure 27. A comparison of the power spectrum of density perturbations in a 
global texture model (solid line) and in the standard CDM model (dashed curve). 
Note the larger power in the texture model on large scales (from Ref. 143). The 
units of P(k) are arbitrary. 

As in the cosmic string model, also in the global texture scenario the length 
scale of the dominant structures is the comoving Hubble radius at t eq . Textures 
generated at t eq are the most numerous, and the perturbations induced by them 
have the most time to grow. 

In both cosmic string and texture models, the fluctuations are non-Gaussian, 
which means that the Fourier modes of the density perturbation 5p have non- 
random phases. Most inflationary Universe models, in contrast, predict (in linear 
theory) random phase fluctuations which can be viewed as a superposition of small 
amplitude plane wave perturbations with uncorrelated phases (for some subtle is- 
sues see Refs. 144 and 145). 
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8. Comparison 



8.1. Role of Hot and Cold Dark Matter 

Before discussing some key observations which will allow us to distinguish 
between the different models, I will discuss the role of dark matter. The key issue 
is free streaming. Recall that cold dark matter consists of particles which have 
negligible velocity v at t eq , the time when sub- horizon scale perturbations can 
start growing: 

v(t eq )<^l (CDM). (8.1) 
For hot dark matter, on the other hand: 

v(t eq ) ~ 1 (HDM) . (8.2) 

Due to their large thermal velocities, it is not possible to establish HDM pertur- 
bations at early times on small scales. Fluctuations are erased by free streaming 
on all scales smaller than the free streaming length 

X c j(t) = v(t)z(t)t (8.3) 

(in comoving coordinates). For t > t eq , the free streaming length decreases as 
£-1/3 rj^g max j ma ] streaming length is 

A7 X = We q ) (8.4) 

which for v(t eq ) ~ 0.1 (appropriate for 25 eV neutrinos) exceeds the scale of galax- 
ies. 

In inflationary Universe models and in the texture theory, the density pertur- 
bations are essentially dark matter fluctuations. The above free streaming analysis 
then shows that, if the dark matter is hot, then no perturbations on the scale of 
galaxies will survive independent of larger-scale structures. Hence, these theories 
are acceptable only if most of the dark matter is cold. 
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Cosmic string theories, in contrast, work well - if not even better - with hot 
dark matter 146 ' 147 ' 7 ' 136 ' 137 ) . The cosmic string seeds survive free streaming. The 
growth of perturbations on small scales A is delayed (it starts once A = Aj(t)) but 
not prevented. 

Let us summarize the main characteristics of the cosmic string, global texture, 
and inflationary Universe theories of structure formation. Inflation predicts ran- 
dom phase perturbations. The density peaks will typically be spherical, and the 
model is consistent with basic observations only for CDM. The global texture and 
cosmic string models both give non-random phase perturbations. The topology is 
dominated by spherical peaks for textures, whereas it is planar or filamentary for 
cosmic strings (depending on the small scale structure on the strings). Textures 
require CDM, whereas cosmic strings work better with HDM. 

8.2. Large-Scale Structure Signatures for Cosmic Strings and 
Textures 

The genus curve 148 ) is a statistical measure for the topology of large-scale struc- 
ture. Given a smooth density field p(x), we pick a density po and consider the 
surface S po where p(x) = po and calculate the genus u(S) of this surface 

v = # of holes of S — #of disconnected components of S . (8.5) 

The genus curve is the graph of v as a function of po- 

For perturbations with Poisson statistics, the genus curve can be calculated 
analytically (Fig. 28). The inflationary CDM model in the linear regime falls in 
this category. The genus curve is symmetric about the mean density p. In the 
texture model, the symmetry about p is broken and the genus curve is shifted to 
the left 143 ). In the cosmic string model, there is a pronounced asymmetry between 
p > p and p < p. At small densities, the genus curve measures the (small number) 
of large voids, whereas for p > p the curve picks 149 ) out the large number of high 
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Figure 28. The genus curve of the smoothed mass density field in a cosmic string 
wake toy model compared to the symmetric curve which results in the case of a model 
with a random distribution of mass points. The vertical axis is the genus (with genus 
zero at the height of the "x"), the horizontal axis is a measure of density ("0" denotes 
average density). 

density peaks which result as a consequence of the fragmentation of the wakes (Fig. 
28). 

The counts in cell statistic 150 ) can be successfully applied to distinguish be- 
tween distributions of galaxies with the same power spectrum but with different 
phases. The statistic is obtained by dividing the sample volume into equal size 
cells, counting the number f(n) of cells containing n galaxies, and plotting f(n) 
as a function of n. 

We 151 ) have applied this statistic to a set of toy models of large-scale struc- 
ture. In each case, the sample volume was (150 Mpc) 3 , the cell size (3.75 Mpc) 3 , 
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and the samples contained 90,000 galaxies. We compared a texture model (galax- 
ies distributed in spherical clumps separated by 30 Mpc with a Gaussian radial 
density field of width 9 Mpc), a cosmic string model dominated by filaments (all 
galaxies randomly distributed in filaments of dimensions (60 x 4 x 4) Mpc 3 with 
mean separation 30 Mpc, a cosmic string wake model (same separation and wake 
dimensions (40 x 40 x 2) Mpc 3 ), a cold dark matter model (obtained by Fourier 
transforming the CDM power spectrum and assigning random phases) , and a Pois- 
son distribution of galaxies. 

As shown in Fig. 29, the predicted curves differ significantly, demonstrating 
that this statistic is an excellent one at distinguishing different theories with the 
same power spectrum. The counts in cell statistic can also be applied to effectively 
two dimensional surveys such as single slices of the CFA redshift survey 9 ). The 
predictions of our theoretical toy models are shown in Fig. 30. 

A third statistic which has proved useful 152 ) is distinguishing models with 
Gaussian and non-Gaussian phases is the void probability function p(R), the prob- 
ability that a sphere of radius R contains no galaxies. 

8.3. Signatures in the Microwave Background 

Inflationary Universe models predict essentially random phase fluctuations in 
the microwave background with a scale invariant spectrum (n — 1). Small devia- 
tions from scale invariance are model dependent and have recently been discussed 
in detail in Refs. 105. In all models, the amplitude must be consistent with struc- 
ture formation. COBE discovery 11 ) of anisotropies in the CMB has provided severe 
constraints on inflationary models. They are only consistent with the present data 
if the bias parameter b is about 1, which must be compared to the value b ~ 2.5 
which is the best value for galaxy formation in this model 153 ). Note that full sky 
coverage is not essential for testing inflationary models since in any set of local ob- 
servations of ST/T, the results will form a Gaussian distribution about the r.m.s. 
value. Mixed dark matter models do slightly better and have recently been studied 
vigorously. 
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Figure 29. The three dimensional counts in cell statistics for a Poisson model (G), 
a cold dark matter model (CDM), cosmic string wakes (SW), string filaments (SF) 
and textures (T). 

Cosmic string models predict non-Gaussian temperature anisotropies. One 
mechanism gives rise to localized linear temperature discontinuities 154 ) ; its origin 
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Figure 30. The two dimensional counts in cell statistic for a slice of the Universe of 
the dimensions of a CFA slice, evaluated for the same models as in Fig. 32. 

is illustrated in Fig. 31. Photons passing on different sides of a long straight string 
moving with velocity v reach the observer with a Doppler shift 
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Figure 31. Sketch of the mechanism producing linear discontinuities in the mi- 
crowave temperature for photons 7 passing on different sides of a moving string S 
(velocity v). O is the observer. Space perpendicular to the string is conical (deficit 
angle a). 

5T 

— ~ 8irG[xv7(v) . (8.6) 

To detect such discontinuities, an appropriate survey strategy ( e.g., full sky survey) 
with small angular resolution is crucial. The distribution of strings also gives rise 
to Sachs- Wolfe type anisotropics 155 ) . 

The theoretical error bars in the normalization of CMB anisotropics from 
strings are rather large - a direct consequence of the fact that the precise form 
of the scaling solution for the string network is not well determined. Nevertheless, 
we can consider a fixed set of cosmic string parameters and ask whether the nor- 
malizations of Gfx from large-scale structure data and from COBE are consistent. 
This has been done numerically in Ref. 106, and using an analytical toy model in 
Ref. 108. 

The analytical model 108 -* is based on adding up as a random walk the individual 
Doppler shifts from strings which the microwave photons separated by angular scale 
v pass on different sides, and using this method to compute AT/T(9). Using the 
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Bennett-Bouchet 156 ) string parameters, the result for G\i becomes 

G\t = (1.3±0.5)l(r 6 , (8.7) 

in good agreement with the requirements from large-scale structure formation 7 ). 

To detect the predicted anisotropics from textures, it is again essential to have 
a full sky survey. However, angular resolution is adequate this time, since the 
specific signature for textures is a small number (~ 10) of hot and cold disks with 
amplitude 157 -* 

<Trp 

— ~ 0.06 x 16 7T Grf ~ 3 • 1(T 5 (8.8) 

and angular size of about 10°. The hot and cold spots are due to photons falling into 
the expanding Goldstone boson radiation field which results after texture collapse 
or due to photons climbing out of the potential well of the collapsing texture 158 ) 
(see Fig. 32). 



Figure 32. Space-time diagram of a collapsing texture (backward light cone) and 
the resulting expanding Goldstone boson radiation (forward light cone). Unwinding 
of the texture occurs at point "TX" . The light ray 72 falls into the potential well and 
is blueshifted, the ray 71 is redshifted. 

Note that the texture model is not ruled out by the recent COBE results. The 
amplitude (8.8) is lower than the pixel sensitivity of the COBE maps. However, the 
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predicted quadrupole CMB anisotropy (normalizing 77 by the large-scale structure 
data) exceeds the COBE data point by a factor of between 2.5 and 4 107 ' 1 59) _ jj enc6) 
biasing must be invoked in order to try to explain the large-scale structure data 
given the reduced value of rj mandated by the discovery of CMB anisotropies. 

8.4. Conclusions 

Topological defect models of structure formation generically give rise to a scale 
invariant power spectrum and are hence in good agreement with the recent COBE 
results on anisotropies of the CMB. The amplitude of the quadrupole temperature 
fluctuation can be used to normalize the models. For cosmic strings, the resulting 
normalization agrees well with the normalization from large-scale structure data. 
For textures, there is a mismatch which requires introducing biasing. For textures, 
the situation is comparable to that in the CDM model, where COBE demands a 
bias parameter b ~ 1, whereas galaxy formation is said to demand 153 -* b ~ 2.5. 

It was emphasized that r.m.s. data intrinsically is unable to differentiate be- 
tween topological defect models (with non-random phases) and inflationary modes 
(with random phases). We need statistics which are sensitive to nonrandom phases. 

The most economical model for structure formation may be the model based 
on cosmic strings and hot dark matter. It requires no new particles (although it 
does require a finite neutrino mass), it agrees well with COBE and with the CFA 
redshift data, and it has clear signatures both for large-scale structure and CMB 
statistics. 
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9. Microphysics of topological defects 



In Chapters 5-8, applications of topological defects were considered for which 
only the gravitational effects of such defects were important. However, there is 
much more structure in defects than is visible through gravity only. Topologi- 
cal defects are localized coherent excitations of gauge and Higgs fields. Matter 
interactions with the defects can probe this structure, can give rise to effects of 
potential importance in cosmology, and possibly even to new direct signatures for 
defect models. In this and the following chapter, I will briefly touch on two aspects 
involving microphysics. The discussion is not intended to be complete, but rather 
wet the curiosity of the reader. First, we consider scattering of fermions by strings 
and monopoles. 

9.1. Callan-Rubakov Type Effects 

Consider cosmic strings arising in a GUT model. The gauge and Higgs fields 
excited in the string mediate baryon and lepton number violating interactions, e.g. 
(in the case of the Higgs field 0) via a term Cj in the interaction Lagrangian 

= (9- 1 ) 

where g is the Yukawa coupling constant and ip is a G multiplet of fermion fields, 
G being the gauge group. 

Due to the term (9.1) in the interaction Lagrangian (and similar terms due to 
the GUT gauge fields coupling to t/j) there is a nonvanishing cross section a for a 
quark to be scattered into a lepton by the defect. The conserved charges which 
are different in the initial and final state are absorbed by the defect (which is here 
treated as a static background field). 

In the case of fermion-monopole scattering, we would naively expect the cross 
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section to be the geometrical one a, 



Vgeom ~ V 2 , (9-2) 

i]^ 1 bein the monopole radius. However, as first discovered by Callan 16 °) and 
Rubakov 161 ) , there is in fact an enhancement of the cross section and 

a~m~ 2 , (9.3) 

m being the fermion mass. Note that (9.3) represents a strong interaction cross 
section. Is there a similar enhancement for cosmic strings? 

Classical physics considerations indicate that there will be no enhancement 
of the baryon number violating inelastic cross section for fermion-cosmic string 
scattering 162 ). For monopoles, the enhancement of the cross section can be viewed 
162 ) as the consequence of an attractive magnetic moment-magnetic field force for 
s wave scattering (see Fig. 33). The potential energy is 

V(r)~-fa.B)~^, (9.4) 

where B_ is the U(l) magnetic field of monopole, and \x is the magnetic moment of 
the fermion. This results in an attractive force 

F(r) ~ 1 . (9.5) 

In the case of ordinary (i.e. non-superconducting) cosmic strings, there is no long 
range magnetic field and hence no enhancement of the inelastic scattering cross 
section due to classical physics effects 162 - ) . For superconducting cosmic strings, B_ 
is circular in the plane perpendicular to the string, and hence again V(r) vanishes 
for s wave scattering 163 \ 
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Figure 33. The magnetic field B of a magnetic monopole M is parallel to the 
magnetic moment \x of an incident s-wave fermion /. 

However, this is not the full story. Quantum mechanical effects due to the non- 
vanishing gauge field at long distances from the string may induce an enhancement 
of the cross section, an effect reminiscent but not identical to the Aharonov-Bohm 
164 ) effect. (The AB effect appears in elastic scattering whereas our effect is due to 
inelastic scattering). Let us therefore consider the quantum mechanical scattering 
of fermions by cosmic strings. 

Fermion-cosmic string scattering was first considered using a first quantized ap- 
proach in Refs. 165 and 166; and extended to inelastic scattering in Ref. 167. The 
same results can be derived 168 ) using a perturbative second quantized approach 
discussed in Ref. 162. 

We consider fermion-defect scattering to first order in perturbation theory, 
treating the defect fields as a fixed background. The transition matrix element is 



where D stands for the defect and S is the scattering matrix. To first order in g 



A = oo< 4> f D\tpD > 



—oo — 



< i)'D\S\il)D > 



(9.6) 





(9.7) 
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where we have split the Hilbert space expectation value into its two tensor com- 
ponents. 

From A we obtain the differential cross section by the usual relation 

% ~ h v j ^ • < 9 - 8 > 

where k' is the momentum of the final particle and V, T are cutoff volume and time 
respectively. 

The procedure 162 ) now is to first evaluate A using free fermion wave functions, 
obtaining the unamplified cross section (da/dVt) . In a second step, we solve the 
Dirac equation in the presence of the defect and determine the amplification factor 

A = (99) 

IpfreeiV *) 

as the ratio of the actual and free wave functions tp and ipf ree , respectively, eval- 
uated at the core radius. Care must be taken to normalize both wave functions 
in the same way at infinity. Note that r = (or p = 0) is taken to be the defect 
center for monopoles (or strings). The actual cross section becomes 

da 4 f da\ 



(four powers of A because in \A\ 2 two fermion wave functions appear and both are 
squared) . 

As an example 168 ) , consider the defect to be a cosmic string. The "geometric" 
amplitude Ao for a string along the z axis is 

A ~ nnn~ 2 dtdz — i{E k -E k ,)t i{k' z -k z )z /q -i \ 

Ao grjv J Mdz {EkEki)1/2 yi/2 e e V- n ) 

The first factor of rj comes from the amplitude of <fi inside the core, integration over 
the core gives r]~ 2 , the factor m comes from fermion spin sums (strictly speaking the 
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summation only occurs when computing the cross section), and the other factor 
m/i] is due to the suppression of the free s-wave amplitude at the core radius. 
Taking the square and performing one z and t integral gives 

(A) 2 ~ 9 2 (^j 2 ^y^(E k - E k ,)8{k z - k z ,) (^j 2 (9.12) 

The two delta functions express conservation of fermion energy and momentum 
along the string. L is a cutoff length. Integrating over k!_ and taking Eft ~ m gives 
the following differential cross section per unit length 

da o f m\ A 1 



*ML 9 UJ m" (9 ' 13) 
The second step of the analysis is to calculate the amplification of the fermion wave 
function at the core radius by solving the Dirac equation in the background field 
of the string: 

ipip-m^ = (9.14) 

The covariant derivative is 

D ll = d ll + ieA ll . (9.15) 

Here, e is the charge of the fermions with respect to the subgroup of G for which 
is the gauge field. We assume that in units of 27r/e the flux of the string is a. 

The solution of the Dirac equation makes use of the techniques pioneered by 
de Vega 169 ). Since we are interested in field configurations independent of z, we 
choose a basis of 7 matrices in which only j z mixes the upper and lower 2-spinors 
of the 4-spinor ip. 

*-(£) (9 - l6) 

(see Refs. 162 and 168 for details). The Dirac equation (9.14) then decouples into 
two separate equations for ip u and ipd respectively. Each is a system of two coupled 
first order differential equations. 
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Next, we expand ip u (and similarly ip^) in partial waves 



n=— oo 



oc 



( 



G n (p)e ; 



) 



irvd—iut 



(9.17) 



where p and i? are the polar coordinates in the plane perpendicular to the string. 
The radial function F n (p) satisfies a Bessel equation with index depending on a. 
It is necessary to solve the equation inside and outside the core, and to match the 
solutions at the core radius. For simplicity, we take the flux of the string to be 
concentrated on the cylinder p = p^ 1 . 

Finally, for given a we search for the mode n(a) for which F n (r7 _1 ) is largest 
and use this coefficient to determine the amplification factor A. We find nontrivial 
amplification: 



where p(a) is plotted in Fig. 34. 



Figure 34. The amplification exponent p plotted as a function of the fractional flux 
a on the string, in a model in which the fermions do not couple to the Higgs field. 

We conclude that in general, there is enhancement of the inelastic baryon 
number violating cross section due to intrinsically quantum mechanical effects, 




(9.18) 
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namely the presence of a nonvanishing gauge potential in the far field limit. The 
maximal cross section is 



da 



dQdL 



~ — , (9.19) 
m 



i.e. a strong interaction cross section. However in general da/dVLdL is suppressed, 
the suppression rate depending on the flux/charge ratio a. 

9.2. NONABELIAN AHARONOV-BOHM EFFECT 

When fermions are scattered by a flux tube, the elastic cross section is also 
much larger than the geometric cross section, as was first discussed by Aharonov 
and Bohm 164 \ The scattered wave has amplitude 

•d being the scattering angle, and k the wave number of the fermion. 

The above corresponds to scattering of a fermion by a U(l) cosmic string. It is 
not hard to see what happens when extending this analysis to fermions scattering 
by a GUT string 170 ). 

Consider a GUT theory with a nonabelian string, i.e. a string in which the 
generator M of the gauge field excited in the string is nondiagonal in the basis 
in which quarks and leptons are separate entries in the fermion multiplet. As an 
example 171 ^ 174 ) ) consider G = SO (10) with the first stage of symmetry breaking 

50(10) -> SU(5) x Z 2 (9.21) 

giving rise to strings. The fermions are in the 16 representation of 50(10) . We 
choose a basis in which 

^ = (ui, U2, «3, "e, d\, d 2 , d 3 , e", df , d c 2 , d% e + , -u\, -u c 2 , -u% -i%) (9.22) 

where U and d denote the usual quarks and the subscripts 1 — 3 stand for the three 
colors. 
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A string solution exists 174 ) for which 
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(9.23) 



This is a nonabelian string for which mixes quarks and leptons. We wish to 
consider the elastic scattering of fermions for this string 17 °). 



We may focus on a two-dimensional subspace, e.g. the us, v e subspace for 
which 

'0 1 
1 



M = 



(9.24) 



To reduce the scattering problem to that for nonabelian string we diagonalize M. 
In the diagonal basis, the fermion eigenstates are 

^1 = ^= ( M 3 + v e ) 



^ 2 = —( U3 -is e ). (9.25) 

The charges are opposite. Hence, the scattering amplitudes differ. If the incident 
beam of fermions is pure M3, the scattered beam will be a mixture of M3 and u e . 
Hence, we have a nonabelian Aharonov-Bohm (AB) effect. The scattering violates 
baryon number! 

The above discussion (see also Refs 175 and 176 for an SU(2) example) has 
been sloppy. Account must be taken of the fact that the physical basis (i.e. the 
basis in which a given entry in ip has fixed physical charges) rotates as we go 
around the string. A recent careful analysis by Ma 177 ) shows that baryon number 
violation in the nonabelian AB effect persists. 
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10. Topological Defects and Baryogenesis 



Even through baryon number violation processes can be mediated by scattering 
of fermions from topological defects, no net baryon asymmetry can be created 
because both initial and final particles are in thermal equilibrium. As realized 
already by Sakharov 178 ), in order to generate a nonvanishing baryon to entropy 
ratio tib/ s, it is necessary to have - in addition to the presence of baryon number 
violation processes - CP violation, and the processes which violate baryon number 
must occur out of thermal equilibrium. Two ways in which topological defects can 
contribute to baryogenesis are sketched in this chapter. The first recurs at the 
unification scale, the second at the electroweak scale. 

10.1. Collapsing Cosmic Strings and Baryogenesis 

The standard mechanism 179 ) of GUT scale baryogenesis is based on the out- 
of-equilibrium decay of superheavy Higgs and gauge particles which freeze out of 
the initial thermal bath of particles at a temperature Tf of the order of the critical 
temperature T c of the symmetry breaking phase transition. 

The above is a viable mechanism. Baryon number violating interactions arise in 
the Lagrangian (see (9.1)), CP violation already exists in the standard model, and 
the out of thermal equilibrium condition is achieved as described in the previous 
paragraph. 

However, if Tf is significantly lower than the mass mx of the superheavy fields, 
then the net baryon asymmetry from the above process may be exponentially 
suppressed since 18 °) for Tf < T < mx the number density nx of X particles 
decays exponentially 

nx ~ e t~ . (10.1) 

As first pointed out in Refs. 181 and 182, cosmic strings can contribute to baryoge- 
nesis in the subset of GUT models which admit cosmic string solutions. The topo- 
logical defects form out of equilibrium concentrations of GUT Higgs and gauge 
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fields. Any emission of quanta form the strings proceeds through ub violating 
interactions, and hence a net baryon number may be generated. In Ref. 181, 
baryogenesis from cosmic strings collapsing onto a line was considered, and in Ref. 
182 the contribution of cusp annihilation 183 ) to ns/s was calculated. A third 
source 184 ) is the final collapse of a cosmic string loop. 

As discussed in Chapter 4, cosmic string loops slowly shrink by emitting grav- 
itational radiation. Once the radius R becomes comparable to the thickness w 
there is no topological barrier for the remnant ball of energy 

E ~ 2itnw ~ 2ttA~ 1/2 ?7 (10.2) 

to decay into Nq quanta of superheavy particles which in turn decay into low mass 
particles producing a net rig. From (10.2) and mx ~ rj it follows that 

N Q ~ 2ttA _1 . (10.3) 

If AB is the average net baryon number per decay, then the net baryon asymmetry 
produced per loop is 

AN B ~ 27rA _1 AB . (10.4) 

Integrating over the distribution of loops after the time of the phase transition and 
dividing by the entropy density we obtain 184 ) 

— ~ X 3 N~ 1 N Q AB (10.5) 

S 

where A is the self coupling constant of the Higgs field potential and iV is the 
number of spin degrees of freedom in thermal equilibrium. In order to obtain this 
result we took Tf = T c and the mean separation of strings to be given by (4.9). 

Since the number density of string loops decays only as a power law (redshift 
effect), the contribution of collapsing cosmic string loops to ns/s decays only as 
Tp 3 and hence may dominate if Tp mx- 
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The result (10.5) must be compared to the maximal ub/s which can be pro- 
duced by the usual GUT mechanism 185 ); 

— ~ N^NxAB . (10.6) 
s 

where Nx is the number of helicity states of all superheavy particles contributing 
to AB. Note that AB depends on the particular GUT model chosen. It enters 
both the standard and topological defect baryogenesis mechanism in the same way. 
The value of AB has been estimated 185 - ) to lie in the range 10~ 13 < AB < 10" 2 . 
Thus, it is in principle possible to generate the observed ub/s ratio which lies in 
the range 

^ ~ 10- 10 . (10.7) 



10.2. Electro weak Baryogenesis with Electro weak Strings 

Recently, it has been realized that nonperturbative processes associated with 
the SU (2) anomaly 186 ) of the standard model can erase a primordial B + L asym- 
metry above the critical temperature of the electroweak symmetry breaking tran- 
sition. Thus, unless a primordial B — L asymmetry is generated at the GUT scale 
(or extra symmetries such as those arising in supersymmetry are introduced 187 - ) ), 
the observed ns/s must be produced below the electroweak scale. 

Initial ideas on how to generate a nonvanishing ns/s were discussed in Ref. 
188, and recently several concrete mechanisms were suggested 189 ~ 191 ). All of 
these mechanisms only work if the electroweak phase transition is first order and 
therefore proceeds by bubble nucleation. 

In electroweak baryogenesis, ng violation recurs nonperturbatively through 
the anomaly, CP violation is already contained in the standard model (although 
stronger CP violation is introduced in Refs 189-191 by enlarging the Higgs sector 
of the model). The role of the first order phase transition in Refs. 189-191 is to 
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generate bubble walls in which the ns violating processes can take place out of 
thermal equilibrium. 

However, the order of the electroweak phase transition is not yet known. There- 
fore it is important to point out that electroweak baryogenesis mechanisms also 
exist if the phase transition is second order 192 ) . Topological defects may play the 
role of the bubble walls in a first order transition. 

However, in the electroweak model and minimal extensions thereof such as 
the two Higgs doublet model there are no stable topological defects. But, as 
realized recently by Vachaspati 193 - ) and considered much earlier by Nambu 1M \ 
energetically stable string-like defects can be constructed in the standard model by 
embedding the Nielsen-Olesen U(l) string into the standard model in the following 
way: 

* = fso(p)j* (^J (10.8) 

Z = A NO (10.9) 

W' = W 2 = A = , (10.10) 

where /noCp)^ and Ano are the Nielsen-Olesen string configurations. 

It can be shown 195 ) that for large Weinberg angles, the above electroweak 
string configuration is energetically stable. Generally, electroweak strings have 
finite length 193 ). Their ends can be viewed as a monopole-antimonopole pair. 

Electroweak baryogenesis works as follows 189 > 190 ). In the outer shell of the 
bubble wall where 

m<gv, (io.il) 

g being the gauge coupling constant, sphaleron 19 °) and local texture unwinding 
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189) even t s take place out of thermal equilibrium at a rate 



T B ~ a^T 4 . (10.12) 

Since there is CP violation in this region, a net baryon asymmetry will be generated. 
It is important that the region of n B violation is in motion with a preferred direction 
(otherwise there is no net CP violation). 

In a model with a second order phase transition and electroweak strings, the 
role of the bubble wall is played by the tips of the string which are moving towards 
each other due to the tension of the string 192 ) . 

The typical length and separation of strings immediately after the phase tran- 
sition is the Ginsburg length (see (4.9)) 

£(t G ) ~ A-V 1 • (10.13) 

The net baryon asymmetry produced by one electroweak string is 

ta+Ats 

/IT 7 
dt—At c eT B , (10.14) 

to 

where dV/dt is the rate of change of the volume in which n B violating processes 
are taking place out of equilibrium, At c is the length of time any point in space 
remains in the above volume, e is the CP violation parameter defined such that 
sYb is the net rate of baryon number producing processes per unit volume. Ats 
is the contraction time of the string. 

Using (10.12) and the geometrical values for dV/dt,At c and At$ (see Ref. 192 
for details), we obtain 

^ ~ \g 2 ea w A , (10.15) 

which is smaller than what is obtained with a first order phase transition only be a 
geometrical factor (the region in which ns violation takes place out of equilibrium 
is much smaller). 
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Note, however, that the above model of electroweak baryogenesis using topo- 
logical defects arising in a second order phase transition only works if extreme con- 
ditions on the parameters of the theory are satisfied: we require stable electroweak 
strings, a second order transition, sufficient CP violation (e ~ 1), constraints on 
the coupling constants to make sure that sphalerons or textures fit into the re- 
gion |0| < gi] and that violating effects in the true vacuum are suppressed at 
to + At s . Nevertheless, it is an existence proof that electroweak baryogenesis with 
a second order phase transition is possible. 
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